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1.0 Introduction 


1.1 Statement of Work 

The work involved in the current task is a continuation of work under a 
previous contract with MSFC. In the first contract a tether simulation was 
obtained from Analytical and Computational Mathematics (ACM) of Houston, Texas. 
It was implemented on the MSFC Sigma V computer. At that time several errors 
were found. These were corrected and documented in the final report for that 
contract. 


The current effort has been extended to include more tasks, but it is 
the purpose here to discuss only the original statement of work for the con- 
tinuation of tether analysis. 

Task A. For a range of tether lengths, end masses, and orbits define/analyze 

tether deployment concepts from the Orbiter (tether motion, payout rate, length, 
tension, end mass position and disturbances) for steady state/dynamic and up/ 

down deployments and from ci rcular/el i ptical orbits. 

Task B. For a range of tether lengths, end masses, and orbits define/analyze 

end mass releasing concepts (imparted delta velocity and final orbits) for 
steady state and dynamic releases taking into account tether and end mass motion 
before and after release. 

Task C. For a range of tether lengths, end masses, and orbits define/analyze 

tether retrieving or disposing concepts (tether motion, reel-in rate, length, 
and Orbiter position) for both reusable and disposable tethers. 

Task D. Perform special tether analysis tasks. 

Task E. Instal 1 /update existing tether programs on the MSFC VAX 11/780 com- 
puter. 

Task F. Submit monthly letter reports and a final report and give a midterm and 
final review. 


1.2 Overview of Work 

The ACM tether simulation was brought up to speed on the HP 9000 .located 
at Control Dynamics. During the verification process, errors were discovered in 
the code and corrected. This version was then implemented on the PD VAX 11/780 
at MSFC as Version 1.0. A few more code errors were discovered and also cor- 
rected. A 4th order Runge-Kutta integration scheme and an initial attempt at 
constant tension deployment were added to the simulation and delivered as 
Version 2.0 at MSFC. MSFC implemented formatting changes for the output files 
and this became Version 2.1. A more refined tension controlled deployment 
scheme was developed, as were terms for material damping in the tether. These 
changes comprise Version 3.0 currently existing at MSFC. 
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A parallel effort to the work above was to update the simulation plot- 
ting capabilities on the MSFC Tektronics graphics terminals. A menu type program 
was developed from which the user can specify types of plots desired, plot 
limits, graph titles, and scale factors. Walking plot graphs were also devel- 
oped which depict the tether motion as it deploys. 

Control Dynamics gave the midterm presentation on February 7, 1986. 
This presentation included the work performed on the simulation through Version 
2.1 and the status of the graphics. 

A complete listing of the code would violate the license agreement 
between ACM and Control Dynamics, so only those lines originating with or modi- 
fied by Control Dynamics can be published. Contact Control Dynamics or the 
Marshall Space Flight Center, Attention C. Rupp PS04, MSFC, AL 35812, regarding 
the license agreement. 
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2.0 Simulation Changes 


2.1 Corrections to Simulation 

As discussed in the overview, several errors in the original simulation 
were discovered and corrected. Appendix A shows the original lines of code, the 
corrections made, and the reasoning behind the changes. 

These corrections involved changing integer loop variables which had 
been previously defined, removal of redundant steps, and corrections of dimen- 
sions and units. 


2.2 Simulation Deletions 

For Version 3.0 delivered at MSFC it was decided to delete all lines of 
code which had been commented out and never used during the first year of this 
project. Appendix B lists those lines which have been deleted and the subrou- 
tines in which they had appeared. 


2.3 Simulation Additions 

During the past year there have also been numerous miscellaneous addi- 
tions to the simulation not covered under the corrections to the simulation or 
under the topics of tether tension deployment, damping or the plotting addi- 
tions. These additions are listed in Appendix C. The main addition involved 
incorporating a second integration scheme. A 4th Order Runge-Kutta integration 
routine was incorporated to help in decreasing run times or when it is desired 
to have a constant step size. 

An addition which was requested by MSFC was to have several orbital 
parameters included in the LGIBLE file output. This change was made in the 
Subroutine Uprint and the variables are the semi -major axis, the eccentricity, 
the argument of perigee, the inclination, the ascending node, the true anomaly, 
the perigee, the apogee, and the longitude of perigee. These terms are all 
derived from the EVECEL(6) vector calculated in Subroutine Vecele. 

Capabilities were also added to Subroutine Deplex to tell the user if 
the tether is fully deployed or fully retracted. If the tether length has 
reached either of these limits, then the deployment scheme is switched to con- 
stant length deployment and the deployment velocity and acceleration are set to 
zero. 

There were also numerous changes made to the format statements in 
Subroutine Uprint. MSFC made these changes internally and relayed them to 

Control Dynamics to implement in the Control Dynamics' subsequent versions of 
the simulation. 
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2.4 Revised Subroutine Descriptions 


Tether (Main) main driver, checks for restart case 


Initialization 

Ubegin 

Uiniti 

Rbegin 

Riniti 

Coot 

Earth 

Uxtoa 

Rotsys 


initializes flags and constants; thrust maneuvers if required 

initializes integration variables 

same as Ubegin for restart case 

same as Uiniti for restart case 

sets mathematical and physical constants 

sets zonal and harmonic terms for non-spherical earth 

rearranges integration variables for integrator 

transformation matrix [RI] (inertial to rotating) 


Program Control 
Ugen 
Timing 
Switch 
Ustop 
Trackf 
Uprint 
Vecele 
Rwrite 


control for integration loop 
time control for thrust maneuvers 
thrust maneuver control 
sets stop flag 

follows max and min tension on tether 

routine to print data to output, plot and restart files 

computes true anomaly for output 

prints data for restart input 


Integration Loop 

Uint driver 

RK78 Runge-Kutta-Fehl berg 7(8) routine and 4th order Runge-Kutta 

routine 


Derivative Calculations 
Uder driver 

Uatox places variables in form used to calculate forces 

Help sets intermediate values for derivative evaluation 

Deplex computes time dependent quantities according to deployment laws 
Tendpl computes tether velocity and components of tether acceleration 

for tension deployment 


Force calculations: 


Ki nema 
Centra 
Bendi n 
Gi nteg . 
Ghalfi 
Thrust 
Normal 
I next 
Homoge 
Fulext 
Basmat 
Li nequ 
Pertur 
Unosph 


kinematic terms 
central force field 
bending stiffness 

intermediate values used in Bendi n 

II 

forces due to thrust maneuvers 
tether tension 

inextensible tether 
homogeneous tether 
fully extensible tether 

called by Inext or Homoge to create intermediate matrix 
called by Inext or Homoge to solve intermediate matrix 
external perturbations 

perturbations due to non-sphericity of earth 
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Erdpot 



Solun 
Kepequ 
Ngrav 
Vites 
Ai rden 
Dens it 
Sol con 
Eclips 
Rot 

Radi at 



standard routine for calculation of higher order 

zonal and harmonic terms 

third body (sun and moon) perturbations 

intermediate routine to prevent calculating positions 

more than once every time step 

analytical ephemeris of sun and moon 

solves Kepler equation 

non-gravitational terms (drag and radiation forces) 

velocity with respect to earth's atmosphere 

air density from Jaccia model 

air density from simple exponential model 

momentum flux due to solar radiation 

sets eclipse flag 

transformation matrix [TI] (inertial to tether) 
radiation pressure 
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2.5 Revised Variables List 


Variable 

Common 

Variable 

Common 

* A1 

Atmos 

Ipyes 

Dripri 


Normco 

I rad 

Force 

* A2 

Atmos 

Iresta 

Restar 


Normco 

Irsavl 

Scount 

A3 

Atmos 

Irsav2 

Scount 

Alpha 

Depl aw 

Islack 

SI ack 

* B1 

Atmos 

Isun 

Thi rdb 


Normco 

Iteg 

Graton 

* B2 

Atmos 

Ithi rd 

Force 


Normco 

lyes 

Dripri 

* B3 

Atmos 

Kstep 

Restar 


Normco 

Latt 

Atmos 

Beta 

Tespec 

Lone! a 

Tespec 

Bstiff 

Bend 

Nactiv(2) 

Thrus3 

Cl 

Normco 

Ndeq 

Gen 

C2 

Normco 

Nend(2) 

Thrus3 

Dac 

Tdeppa 

Nfmas 

Traqua 

Dacini 

Tindpa 

Nf mi s 

Traqua 

Deg 

Cbasic 

Nodes 

Tindpa 

Direc(2,3) 

Thrus2 

Nprint 

Dripri 

Dt 

Umesh 

Npuls(2) 

Thrusl 

Dvn 

Tdeppa 

Nsat(2) 

Thrus4 

Dvnini 

Ti ndpa 

Nstep 

Gen 

Em 

Atmos 

Nstop 

Stopva 

Epsj2 

Oblate 

Nswitc(2) 

Thrus3 

Fc 

Depl aw 

Ntess 

Oblate 

Formas 

Traqua 

Nthrus 

Thrusl 

Formis 

Traqua 

Ntype{2) 

Thrus2 

Fr 

Depl aw 

Nzon 

Oblate 

Frk 

Depl aw 

OmqSOr 

Cdynae 

Forpro(19) 

Projec 

0mt50r 

Cdynae 

Gamma 

Tespec 

P(19) 

Cocoef 

G1 

Bouwei 

Pa ram 

Gen 

G2 

Bouwei 

Pdurat(2) 

Thrusl 

I bend 

Bend 

Phi 

Atmos 

Idepl 

Depl aw 

Pi 

Cbasic 

Idurn 

Atmos 

Pih 

Cbasic 

Idrag 

Force 

Pi nter(2) 

Thrusl 

Imoon 

Thi rdb 

• R 

Cart 

Inosph 

Force 

Rad 

Cbasic 

Ipl ot 

Dripri 



Iplt 

Dripri 




* Variables that occur in two common blocks never occur in the same sub- 
routine. They represent two different values. In this sense, they act 
as dummy variables. 
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Variable 

Relvel(19,3) 

Relvqu(19) 

Rho 

Rnode(19,3) 

Rnodea(19) 

Rshu 

Rshu2 

Rshutt(3) 

Sl(19) 

S2(19) 

S3(19) 

S4(19) 

Shmass 
Shrefc 
Shucro 
StdSOr 
Step in 
Subcro 
Sumass 
Surefc 
T 

Talpha(2) 

Tau(19,3) 

Taumag(19) 

Tauqua(19) 

Tbeta(2) 

Tdepoc 

Tdnow 

Tedens 

Tess(2,36) 

Tetrad 

Tiipl(19) 

Tini 
Tlaini 
Tlevel (2) 


Common 

Rel CO 

Intern 

Atmos 

Abscoo 

Abscoo 

Shutco 

Shutco 

Shutco 

Cocoef 

Cocoef 

Cocoef 

Cocoef 

Emass 

Radcoe 

Shape- 

Cdynae 

Gen 

Shape 

Emass 

Radcoe 

Umesh 

Thrus2 

Relco 

Interm 

Interm 

Th rus2 

Jdate 

Jdate 

Tespec 

Ckufue 

Shape 

Intern 

Ti ndpa 

Actual 

Thrus2 


Variable 

Tin 

Tlnl 

Tln2 

Tlninl 

Tlnin2 

Tlnini 

Tplast 

Tprint 

Trefc 

Tstart(2) 

Tstop 

Twopi 

Vnode(19,3) 

Vshutt(3) 

X(7) 

Xincre 

XmassO 

Xnassl 

Xmas s 2 

Xmu 

Xmum 

Xnus 

Xnforc(19) 

Y(123) 

Z(123) 

Zl(4) 

Zonal (23) 
Zs(4) 


Common 


Tdeppa 

Retrap 

Retrap 

Retrap 

Retrap 

Tindpa 

Dripri 

Dripri 

Radcoe 

Thrusl 

Stopva 

Cbasic 

Abscoo 

Shutco 

Cart 

Ti ndpa 

Ti ndpa 

Tindpa 

Tdeppa 

Cpoems 

Cpoems 

Cpoems 

Nforce 

Umesh 

Umesh 

Jposms 

Ckufue 

Jposms 
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2.6 REVISED COMMON BLODK LIST 

/ABSCOO/ RNODE(19,3), RNODEA(19) ,VNODE(19,3) 

RNODE(i,j) = component j of the position of the node i 
referred to geocentric inertial reference frame. 

RNODEA(i) = |^i| 

VN0DE(i,j) = component j of the velocity d^^/dt of the node i 
referred to the geocentric inertial reference 
frame. 

These quantities are computed in the subroutine UATOX. 

Used in subroutines: UATOX, UPRINT, CENTRA, PERTUR, NGRAV. 

/ACTUAL/ TLAINI 

TLAINI = deformed length of the tether at initial time. 

Input quantity. 

Used in subroutines: UBEGIN, HEAD, RBEGIN. 

/ATMOS/ PHI, EM, RHO, Al, A2, A3, Bl, B2, B3, IDIURN, LAT 

Constants and flags for the atmospheric model. Input quantities. 
Used in subroutines UINITL, RINITL, UBEGIN, RBEGIN, HEAD, AIRDEN. 

/BEND/ BSTIFF, IBEND 

Flag and constant for the bending stiffness. Input quantities. 

4 

Tia 

BSTIFF=EJ,E=YOUNG's MODULUS J= — 

4 

Used in subroutines: UINITL, RINITL, UBEGIN, RBEGIN, HEAD, UDER, 

BENDIN. 
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/BOUWEI/ Gl, G2 


The weights for the extrapolation formulas of the discretization 
at boundaries. They are set 

Gl = 1.5, G2 = -0.5 

in the subroutines UINITL or RINITL. 

Used in subroutines: UINITL, RINITL, HELP, BENDIN, NORMAL. 

/CART/ X(7), R 


X{l-3) 

= cartesian 

position vector 

X{4-6) 

= cartesian 

velocity vector 

X(7) 

= time 


R 

= magnitude 

of cartesian position vector X(l-3) 


These quantities are set in subroutine PERTUR and used locally 
for the evaluation of the gravitational perturbations. 

Used in subroutines: PERTUR, UNOSPH, UTHIRD. 

/CBASIC/ PI, TWOPI, PIH, DEG, RAD 

Basic mathematical constants. 

PI = IT, TWOPI = 2ir, PIH = ir/2, DEG = 180/ir, 

RAD = ti/180. 

They are set in subroutine COOT. 

Used in subroutines: UINITL, RINITL, COOT, UPRINT, VECELE, 

KEPEQU, THRUST. 



/CDYNAE/ STD50R, 0MT50R, 0MQ50R t 

STD50R, 0MT50R, 0MQ50R = constants to compute the angle W between 

Greenwich meridian and mean equinox as a 
function of the Modified Julian Date in 
radians. 

These quantities are set in subroutine COOT. 

Used in subroutines: COOT, UNOSPH. 

/CKUFUE/ Z0NAL(23), TESS(2,36) 

Z0NAL(l-23) = zonal coefficients of earth geopotential 

TESS(L,l-36) = tesseral coefficients of earth geopotential 
TESS(l,L*L-L)/2+M) = C(L,M) 

TESS(2,(L*L-L)/2+M) = S(L,M) 

These quantities are set in subroutine EARTH. 

Used in subroutines: UINITL, RINITL, EARTH, UNOSOPH, ERDPOT. 

/COCOEF/ Sl(19), S2(19), $3(19), $4(19), P(19) 

Constants and coefficients for system of linear equations. 
Matrix with four non-vanishing diagonals. 

Locally used in subroutines: HOMOGE, INEXT, BA$MAT. 

/CPOEMS/ XMU, XMUM, XMUS 

XMU = gravitational constant of the earth 

XMUM = gravitational constant of the moon 

XMUS = gravitational constant of the sun 

These quantities are set in subroutine COOT. 

Used in subroutines: UINITL, RINITL, COOT, VECLE, CENTRA, UTHIRD. 

t See also: Definition and units in the FORTRAN source listing of the 

subroutine COOT. 
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/DEPLAW/ ALPHA, IDEPL, FC, FRK, FR 


The flag IDEPL specifies the explicit deployment/retraction law 
and if 


IDEPL 
= 1 
= 2 


^(t) = ^(to) = const 

m • 

A(t) = ^(to) = const 


n{t) = ^-(to) = const 


£(t) = a£(t) . (exponential 
.deployment law) 
tension deployment law 


= 3 
= 4 
= 5 

to is the initial time. 

ALPHA = o for exponential deployment law. 
FC = commanded max tension 0 n-1 


FRK = tension coefficient 
FR = tension @ n-1 

Used in subroutines: UINITL, RINITL, UBEGIN, RBEGIN, DEPLEX, 

UXTOA, TENDPL, RWRITE, UATOX, UDER, UGEN, 
UPRINT. 


/DRIPRI/ NPRINT, IPYES, TPLAST, TPRINT, lYES, IPLOT, IPLT 
Output control . 


NPRINT, TPRINT: 
IPYES : 

TPLAST : 

lYES : 

IPLOT : 

IPLT : 


Input quantities. 

flag set >0 in subroutine UGEN if results of the 
just performed integration step have to be printed 

last print time 

stop flag set in USTOP 

counter for plot file 

counter for walking plot file 


Used in subroutines: UBEGIN, RBEGIN, RGEN, UPRINT. 
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/EMASS/SUMASS, SHMASS 


SUMASS 


mass of the subsatellite without retracted part of 
the tether. Input quantity. 

SHMASS 

= 

mass of the shuttle without retracted part of the 
tether. Input quantity. 

Used in 

subroutines: RINITL, UBEGIN, RBEGIN, HEAD, DEPLEX, 

TENDPL. 

/FORCE/INSOPH, 

ITHIRD 

, IDRAG, IRAD 

INOSPH 


determines whether or not non-sphericity pertur- 
bations are to be computed. Input quantity. 

(=1, no; =2, yes) 

THIRD 


determines whether or not third body perturbations 
are to be computed. Input quantity. 

(=1, no; =2, yes) 

I DRAG 

= 

determines whether or not drag perturbations are 
to be computed. Input quantity. 


= 1 

: No perturbations due to drag. 


=2 

: Perturbations using simple air density 

model . 


= 3 

: Perturbations using Jacchia air density 

model . 

I RAD 

= 

determines whether or not radiation perturbations 
are to be computed. Input quantity. 


= 1 

: No radiation perturbations. 


=2 

: Perturbations according to an averaged 

solar. 


= 3 

: Perturbations according to a sinusoidally 

varying solar constant with earth eclipse 
effects included. 


Used in subroutines: 


UINITL, RINITL, UBEGIN, RBEGIN, HEAD, UDER, 
PERTUR, NGRAV. 



/GEN/ STEPIN, NDEQ, 

NSTEP, PARAM 

STEPIN 

= initial guess for the integration step. 
Input quantity. 

NDEQ 

= number of differential equations. Computed in 

subroutine UGEN 

NSTEP 

= counts the number of integration steps. 

PARAM 

= tolerated truncation error on each integrated 
quantity. Input. 

Used in subroutines: UINITL, RINITL, UBEGIN, RBEGIN, HEAD, UGEN, 

UDER, UPRINT, RWRITE, USTOP, UINT. 

/GRATON/ ITEG 

ITEG : 

flag for integration scheme, 1 variable stepsize 
2-4th order Runge - Kutta 

Used in subroutines: RBEGIN, RK78, UBEGIN 

/INTERM/ TAUQUA(19) 

, TAUMAG(19), TIIP1(18), RELVQU(19) 

TAUQUA(i) 

1 1 2 
= lli+1/2 1 

TAUMAG(i) 

= 1 li +1/2 1 

TllPl(i) 

= tli+1/2 • li+3/2] 


RELVQU(i) = I Wi+1^ I ^ 


These quantities are computed in the subroutine HELP and will be 
used for the evaluation of the 2 ,i+l /2 and w^+l^ integration vari- 
abl es . 

Used in subroutines: UPRINT, TRACKF, HELP, BENDIN, THRUST, 

HOMOGE, FULEXT, INEXT, BASMAT, DEPLEX, 

TENDPL, UDER. 


/JDATE/ TDEPOC, TDNOW 

TDEPOC = is the Julian date at initial time (since Jan. 1, 
1950). 

TDNOW = is the Julian date in days after each step of pro- 

pagation. 

Used in subroutines: UBEGIN, RBEGIN, UGEN, RWRITE, USTOP, PERTUR. 
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/JPOSMS/ ZS(4), ZL(4) . 

ZS(l-3), ZS(4) = position and distance of the sun referenced to 
the equatorial system of the earth. 

ZL(l-3), ZL(4) = position and distance of the moon referenced to 
the equatorial system of the earth. 

Used in subroutines: UTHIRD, USOLUN, NGRAV, AIRDEN. 


/NFORCE/ XNF0RC(19) 

XNFORC(i) = N*-j+iy 2 , tether tension magnitude x factor 

Used in subroutines: RBEGIN, UPRINT, TRACKF, RWRITE, NORMAL, 

HOMOGE, FULEXT, INEXT, DPLEX, TENDPL, UDER. 


/NORMCO/ Al, A2, Bl, B2, B3, Cl, C2 

Intermediate coefficients used to calculate normal forces. These 
quantities will be computed in subroutine HELP. 

Used in subroutines: HELP, BENDIN, NORMAL, INEXT, BASMAT. 


/OBLATE/ EPSJ2, NZON, NTESS 

2 

EPSJ2 = 3/2 J 2 r e» computed in subroutine UINITL or 

RINITL. 

where 

J 2 = second zonal coefficient of earth, set in 
subroutine EARTH 

Te = equatorial mean radius of earth, set in 
subroutine EARTH. 

NZON = the desired number of zonal terms to be included 

in non-sphericity model. Input quantity. 

NTESS = the desired number of tesseral terms to be 

included in non-sphericity model. Input quantity. 

Used in subroutines: UINITL, RINITL, UBEGIN, HEAD, UNOSPH, 

RBEGIN. 
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/PROJEC/ F0RPR0(19) 

These are scalar products computed in the subroutine NORMAL and 
used to calculate tether tension. 

Locally used in subroutines: NORMAL, INEXT, BASMAT. 


/RADCOE/ TREFC, SUREFC, SHREFC 

Reflection coefficients for the computation of the radiation 
pressure. Input quantities. 

TREFC = reflection coefficient for the tether 

SUREFC = reflection coefficient for the subsatellite; 

currently not used 

SHREFC = reflection coefficient for the shuttle; currently 

not used 

Used in subroutines: UINITL, RINITL, UBEGIN, RBEGIN, NGRAV. 


/RELCO/ TAU(19,3) RELVEL(19,3) 


These are the integration variables. 


TAU(i,j) = component j of 

RELVEL(i,j) = component j of 

Used in subroutines: UBEGIN, 

RWRITE, 

GINTEG, 

UEPLEX, 


the vector J.i+ly 2 
the vector w^+i/^ 

RBEGIN, UATOX, UXTOA, UDER, UPRINT, 
HELP, KINEMA, CENTRA, BENDIN, 
6HALFI, NGRAV, THRUST, NORMAL, 
TENDPL, KINMAX, FULEXT, HOMOG. 


/RESTAR/ IRESTA, KSTEP 

IRESTA . = determines whether the case to be run is a new 
case or a restart case. Input quantity. 

KSTEP = number of step on the restart file from which a 

restart case has to be initialized. 

If KSTEP does not coincide with a printed integra- 
tion step the program will use the next largest 
printed step for initialization conditions of the 
restart case. 

Used in subroutines: HEAD, RBEGIN. 
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/RETRAP/ TLNINl, TLNIN2, TLNl, TLN2 


TLNINl 

TLNIN2 

TLNl 


TLN2 


undeformed length of that part of the tether which 
is initially retracted in the subsatellite. Input 
quantity. 

undeformed length of that part of the tether which 
is initially retracted in the shuttle. Input 
quantity . 

undeformed length of that part of the tether which 
is retracted in the subsatellite at actual time t. 
Computed in subroutine DEPLEX. Constant in 
current version. 

undeformed length of that part of the tether which 
is retracted in the shuttle at actual time t. 
Computed in subroutine DEPLEX. 


Used in subroutines: UBEGIN, RBEGIN, HEAD, UPRINT, RWRITE 

DEPLEX, USTOP, TENDPL, UDER. 


/SCOUNT/ IRSAVl, 

IRSAV2 





IRSAVl 

= number of rejected 
of the run. 

integration 

steps 

since 

start 

IRSAV2 

= number of rejected 
printout time. 

integration 

steps 

since 

last 


Used in Subroutines: UPRINT, UINT. 


/SHAPE/ TETRAD, SUBCRO, SHUCRO 

The geometrical quantities of this COMMON block will be used for 
the evaluation of air drag and radiation pressure forces. 

TETRAD = Radius of tether cross-section. Input quantity. 

To be given in mm. 

SUBCRO = Cross-sectional area of 2 subsatel 1 ite. Input quan- 

tity. To be given in m . 

SHUCRO = Cross-sectional |irea of shuttle. Input quantity. 

To be given in m . 

Used in subroutines: UINITL, RINITL, UBEGIN, RBEGIN, NGRAV. 
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/SHUTCO/ RSHUTT(3), VSHUTT(3), RSHU2, RSHIJ 

Cartesian coordinates and velocities of the shuttle referred to 
the geocentric inertial reference frame. RSHUTT and VSHUTT are 
integration variables. 

RSHUTT(3) = _rji(t) 

VSHUTT(3) = _^(t) 

RSHU2 = I _T;,(t) I ^ 

RSHU = I rn(t) | 

RSHU2 and RSHU will be computed in the subroutine UATOX. 

Used in subroutines: UBEGIN, RBEGIN, ROTSYS, UATOX, UXTOA, 

UPRINT, RWRITE, KINEMA, CENTRA, PERTUR, 
NGRAV, DEPLEX, TENDPL, KINMAX. 


/SLACK/ ISLACK 

I SLACK is a flag used only in connection with the model for the 
fully extensible tether. ISLACK is an input quantity and if 

ISLACK equal 1 : Negative tether tensions accepted. 

ISLACK not equal 1 : Negative tether tensions not accepted, the 
corresponding normal forces are set zero by the program. 

Used in subroutines: UINITL, RINITL, UBEGIN, RBEGIN, HEAD, 

FULEXT. 


/STOPVA/ NSTOP, TSTOP 

These are input quantities causing a termination of the execution 
when the number of integration steps exceeds NSTOP or when the 
integration time exceeds TSTOP. If negative values are assigned 
to NSTOP or TSTOP, the corresponding stopping condition is not 
active. 


Used in subroutines: UBEGIN, RBEGIN, USTOP 



/THIRDB/ ISUN, IMOON 

Flags to switch on or off the gravitational perturbations due to 
sun and moon and if 

ISUN 

= 0 perturbation due to the sun off 

= 1 perturbation due to the sun on 

IMOON 

= 0 perturbation due to the moon off 

= 1 perturbation due to the moon on 

Used in subroutines: UINITL, RINITL, UBEGIN, RBEGIN, HEAD, 

UTHIRD. 


/THRUSl/ TSTART(2), PDURAT(2), PINTER(2), NTHRUS, NPULS(2) 


Parameters to specify thrust maneuvers. 


NTHRUS 

TSTART(J), 

NPULS(J), 

PDURAT(J), 


: Number of thrust maneuvers. NTHRUS May be equal 

0,1, or 2. Input. 

J=l,2: Start time of thrust maneuvers J. Input. 

J=l,2: Number of pulses of thrust maneuver J. 
Input. 

J=l,2: Duration of one pulse of thrust maneuver J. 

Input. 


PINTER(J), J=l,2: Time interval between subsequent pulses of 

thrust maneuver J. Input. 

Used in subroutines: UINITL, RINITL, UBEGIN, RBEGIN, HEAD, 

TIMING, SWITCH, UDER, THRUST. 
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/TDEPPA/ TLN, DVN, DAC, XMASS2 


Time dependent parameters for deployment or retraction. They 


will be 

computed in 

the subroutine DEPLEX. 

TLN 

II 

Undeformed tetherlength at the time t. 

DVN 

dt 

d^ 

Deployment/ retraction velocity 

DAC 

dt^ 

Deployment/ retraction acceleration 

XMASS2 

= m2(t) : 

Mass of shuttle including retracted part of 
the tether 

Used in 

subroutines: 

UPRINT, RWRITE, DEPLEX, HELP, KINEMA, 
BENDIN, NGRAV, THRUST, NORMAL, HOMOGE, 
FULEXT, INEXT, UXTOA, USTOP, TENDPL, 
KINMAX, RBEGIN, UATOX, UBEGIN, UDER. 

/TESPEC/ TEDENS 

, BETA, GAMMA, LDNELA 

Materi al 

tether quantities specified by input. 

TEDENS 

= u : 

Specific mass of the tether (g/m) 

LONELA 

= Flag 

to choose the model for the longitudinal 


deformation and if LONELA 


= 1 Inextensible tether 


= 2 Homogeneous longitudinal deformation 

= 3 Fully extensible tether. 

BETA = B : Longitudinal stiffness 

= EF E = Young's modulus (Nt/m^) 

2 2 
P = r IT area of tether cross section (m ) 


GAMMA = Damping coefficient 

Used in subroutines: UINITL, RINITL, UBEGIN, RBEGIN, HEAD, UDER, 

UPRINT, DEPLEX, HELP, KINEMA, NGRAV, NORMAL, 
HOMOGE, FULEXT, TENDPL, KINMAX, UATOX, UGEN. 
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/THRUS2/ TLEVEL(2), NTYPE(2), TALPHA(2), TBETA(2), DIREC(2,3) 


Parameters to specify thrust maneuvers. 


TLEVEL(J), J=l,2 : Thrustlevel in Newton. Input. 


NTYPE(J), 

=1 

=2 


J=l,2 : Specifies type of input giving the thrust 
direction and if NTYPE(J) 

: Direction of the thrust force of maneuver J 

given by the direction vector DIREC(J,3). 

: Direction of the thrust force of maneuver J 

given by the angles TALPHA(J) and TBETA(J). 


TALPHA (J), J=l,2: for thrust maneuver J (see Figure 2.6.1). 

Input. 

TBETA (J), J=l,2: B-r for thrust maneuver J (see Figure 2.6.1). 

Input. 


DIREC(J,I), J=l,2, 1=1, 2, 3: Components of the direction of the 

thrust force of maneuver J referred to the 
orbital reference frame. Calculated from 
DIRECl and DIREC2 in UINITL. 


Used in subroutines: UNITE, RINITL, UBE6IN, RBEGIN, HEAD, THRUST. 


/THRUS3/ NACTIV(2), NSWITC(2), NEND(2) 

Flags and counters for thrust maneuver control. 

NACTIV(J), J=l,2: Specifies whether the thruster of maneuver J 

is active or not and if NACTIV(J) 


= 0 


Thruster not active. 

Thruster active. This flag is set in the 
subroutine SWITCH. 


NSWITC(J) J=l,2; 


"switch off" - 
to set the end 
updated in the 


Counts the "switch on" - and 
actions of maneuver J in order 
flag NEND(J). This counter is 
subroutine SWITCH. 

NEND(J) J=l,2: End flag for the thrust maneuver J. The flag 

is initialized to 0 in the subroutines UINITL 
or RINITL and is set equal to 1 in the 
subroutine SWITCH if the thrust maneuver J 
has been completed. 

Used in subroutines: UINITL, RINITL, UBEGIN, RBEGIN, HEAD, UGEN, 

TIMING, SWITCH, UDER, THRUST. 
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(SHUTTLE OB 
SUBSATELLITE) 

Figure 2.6.1. 
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/THRUS4/ NSAT(2) 


NSAT(J) 

J=l,2 : Flag which specifies whether the maneuver J has 
to be applied on the shuttle or on the subsatellite and 
if NSAT(J). 

= 1 

maneuver J to be applied on the shuttle. 

= 2 

maneuver J to be applied on the subsatellite. 

Used in 

subroutines: UINITL, RINITL, UBEGIN, RBEGIN, HEAD, UDER, 

THRUST. 

/TINDPA/ TINI, 

TLNINI, DVNINI, DACINI, XINCRE, XMASSO, 

XMASSl, 

NODES 

TINI 

= to : Initial time for a given explicit deployment/ 
retraction law. In the current version of the 
program. TINI is set zero in the subroutine 
UGEN. (At the same place, the initial value of 
the integration time T is set zero). 

TLNINI 

= T(to): Initial value of the undeformed deployed 
tetherlength. Input. 

DVNINI 

• 

= "5^(to): Initial value of the deployment velocity. 
Input. 

DACINI 

« m 

= X(to): Initial value of the deployment acceleration. 
Input. 

NODES 

= n : The number of spatial nodes including the end- 

masses. Input. 

XINCRE 

= h = l/(n-l) : Increment of the spatial discretization. 

XINCRE is initialized in the subroutine 
UBEGIN or RBEGIN. 

XMASSO 

: Dummy variable. Not used in the current ver- 
sion of the program. 

XMASSl 

= mi(t): Mass of the subsatellite including retracted 
part of the tether. In the current version of 
the program which does not treat deployment 
from the subsatellite, XMASSl is set equal to 
SUMASS In the subroutine UBEGIN. 

Used in 

subroutines: UINITL, RINITL, UBEGIN, RBEGIN, HEAD, UGEN, 

UATOX, UXTOA, UDER, UPRINT, TRACKF, RWRITE, 
DEPLEX, HELP, KINEMA, CENTRA, BENDIN, 

GINTEG, GHALFI, PERTUR, NGRAV, THRUST, 

NORMAL, HOMOGE, FULEXT, INEXT, BASMAT, 

USTOP, TENDPL, KINMAX. 
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/TRAQUA/ FORMAS, FORMIS, NFMAS, NFMIS 


FORMAS : Maximum tethertension with respect to all nodes and 

all integration steps between subsequent print times. 
The quantity is tracked in the subroutine TRACKF. 

NFMAS : Node where FORMAS appears. The quantity is tracked 

in the subroutine TRACKF. 


FORMIS : Minimum tethertension with respect to all nodes and 

all integration steps between subsequent print times. 
The quantity is tracked in the subroutine TRACKF. 

NFMIS : Node where FORMIS appears. The quantity is tracked 

in the subroutine TRACKF. 


Used in subroutines: UPRINT, TRACKF. 


/UMESH/ T, DT, Y(123), F(123) 

Allocation of arrays for the numerical integration. 


T 

DT 

Y(184) 

F(1472) 


Time (independent variable). 

Stepsi ze 

One dimensional array of integration variables. 
Auxiliary array. 


Used in subroutines: UGEN, UPRINT, TRACKF, UINT, DEPLEX, TENDPL. 
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2.7 Revised Input File Example 


Input data for tethered satellites computation 

Iresta Flag to determine input for new or restart case 

F A new run is initialized 

T A restart case is initialized. The state of the 

system is read from the restart file at step>kstep 

Iresta=F,Kstep=4000 


Inosph Flag to switch on/off the perturbation due to the 

nonsphericity of the earth 

1 Off 

2 On 

Nzon Order of expansion of the earth potential in zonal 

terms. Nzon = 0,23 

Ntess Order of expansion of the earth potential in 

tesseral terms. Ntess = 0, ..., 8 

Inosph=l,Nzon=23,Ntess=8 

************************************************************************** 

Ithird Flag to switch on/off the perturbations due to 

third bodies 

1 Off 

2 On 

Isun Flag to switch on/off the perturbation due to the 

sun 

0 Off 

1 On 

Imoon Flag to switch on/off the perturbation due to the 

moon 

0 Off 

1 On 

It hi rd=l, I sun=0, lmoon=0 


Idrag 

1 

2 

3 

Idiurn 

F 

T 


Flag to switch on/off the perturbation due to the 

air drag 

Off 

On, simple exponential density model 

On, sophisticated density model (fit to Jaccia- 

1973) 

Flag to switch on/off diurnal effects for Idrag=3 

Off 

On 
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Latt 


Flag to switch on/off latitudinal effects for 
Idrag=3 
F Off 

T On 

Idrag=3, Idiurn=T,Latt=T 

‘kieificiciTicicieic'kicific:ki(iei(icicieic'kieieie'kieicicic'kic'k-k'kiificic'kicicicieieicicieicicieicic'k'k-k'k'kic'k‘k'kic'k’k'k‘k-k-kic'k'k'k-k 


Tetrad Radius of tether crossection (tnm) 

Shucro Shuttle crossection (ni**2) 

Subcro Subsatellite crossection (m**2) 

Tetrad=.4d0,Shucro=50.d0,Subcro=5.d0 

************************************************************************** 


Phi 

Em 

Rho 

A1,A2,A3, 

B1,B2,B3 


Atmospheric bulge lag angle (deg) 

Model exponent 

Atmospheric density at earth surface (kg/m**3) 
Model fit paramters (fit to Jaccia model) 


Phi=37.d0,Em=2.75d0,Rho=1.225d0 

Al=-.24436605d2,A2=-.14473998d-l,A3=.89553181d3 

Bl=-.21795610d0,B2=.33153996d-2.B3=-.24610244d2 


itificicieiicifiritiric^iciticiticiciiriciricicic^hiciriricirititieici^icicitrieiciirrhiciciiricieieic'hiciciricicifieiciticicicicicic’kiticiririr 


I rad 

1 

2 

3 

Trefc 

Shrefc 

Surefc 


Flag to switch on/off the perturbation due to the 

radiation pressure 

Off 

On, solar constant is averaged 

On, include variations in solar radiation due to 

earth's orbit 

Reflection coefficient for tether 
Reflection coefficient for shuttle 
Reflection coefficient for subsatellite 


Irad=l,Trefc* .8dO,Shrefc= ,8dO,Surefc=.8dO 
************************************************************************** 


Ibend Flag to switch on/off the bending stiffness 

1 Bending stiffness neglected 

2 Bending stiffness included 

Bstiff Rending stiffness (kg*km**3/sec**2) 

Bsti ff=EJ 
E=Young's modulus 
J= Ha**4/4 

a=tether radius (km) 

Ibend=l,Bstiff=5.625d-4 

icicicificicicicicieicicic’kicicicirieiricicicicicicicicicificierkieicificicicicic'k'kieicie-kieieicirit-k'k'k-k’k-k'k'k'k'k-k'kicic'k’kic'k-k-k'kie 
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Lone! a Specifies model for longitudinal deformation 

1 No longitudinal deformation; tether forces 
determined by constraint equations 

2 Homogeneous longitudinal deformation; tether forces 
determined by constraint equations and by Hooke's 
law for the segment next to the shuttle 

3 Extensible tether; tether forces determined by 
Hooke's law 

Beta Longitudinal elasticity (constant in Hooke's law) 

Beta=EA 

E=Young's modulus (Nt/m**2) 

A=cross sectional area of tether (m**2) 

Gamma Damping coefficient 

Islack Slackness flag only for extensible tether 

1 Negative tether tensions accepted 

2 Negative tether tensions not accepted but set to 
zero 

Lonela=l,Beta=5.0d4,Gamma=5.d0, Islack=l 

icieicicicicic'kicicicicicic'kic'kieie-k^ie'kicicic'kicieieiricieicicicicieicic'k-k'k-k'k'k'kic'k'k-kieic'k'k-k'k'k-k-k'k'k'k'k'k-k-k'k-k'kic'k'kic 


Nthrus 
Tstart(j) 
Npuls(j) 
Pdurat( j ) 
Pinter(j) 


Number of thrust maneuvers (0,1,2) 

Start time for thrust maneuver (sec), j=l,2 
Number of pulses for thrust maneuver j 
Duration of pulse for thrust maneuver j (sec) 
Interval between pulses for maneuver j (sec) 


Nthrus=0,Tstart=100.d0,500.d0,Npuls=3,3 
Pdurat=50.d0,50.d0,Pi nter=50.d0,0.d0 


************************************************************************** 


Tlevel (j) 
Ntype(j) 

1 


2 


Nsat(j) 

1 

2 


Thrustlevel (Nt) of maneuver j 
Specifies type of input giving thrust direction 
Thrust direction for maneuvers 1X2 to be given by 
direction vectors Direcl(3) and Direc2(3) 
respecti vely , with respect to the orbital 
reference frame (normalized automatically) 

Thrust direction to be given by angle Talpha(j) 
and Tbeta(j) in degrees with respect to the tether 
direction (see Figure 2.7.1) 

Specifies whether maneuver j is done at the shuttle 

or subsatellite 

Maneuver j at the shuttle 

Maneuver j at the subsatellite 


Tlevel=100.d0,20.d0,Ntype=2,2 
Di recl=0.d0,0.d0,0.d0 
Di rec2=0.d0,0.d0,0.d0 
Talpha=45.dO,0.dO,Tbeta=45.dO,180.dO 
Nsat=2,2 
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z 



FIGURE 2.7.1 


X.Y,Z 



Orbital reference frame 

Projection of subsatellite to YZ-Plane 
Projection of subsatellite to XZ-Plane 
In-plane angle 
Out-of-plane angle 
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Nodes 


Number of spatial nodes including the endmasses 
(4-30) 

Stepin Initial guess for integration time step 

Param Tolerated truncation error for integration 

Nodes=6,Stepi n=l.d-2,Param=l.d-8 


Idepl Specifies tether deployment mode 

1 Constant commanded tether length 

2 Constant commanded deployment/retraction velocity 

3 Constant commanded deployment/retraction 
acceleration 

4 Exponential . law: dVdt = alpha*t 

5 Tension deployment 

Alpha Constant in the exponential law 

Idepl=5,Al pha=l .d-4 

*★**★***★★★★*★**★****•************★*****★*★*★★*•*•**•*★********•**★**★★★★*■★**** 

For tension deployment: 

Absolute value of maximum tension at node next to 
shuttle (Nt) 

Slope coefficient fori tension (Nodes-l)j <Fc 
Initial tension at Nocies-l (Nt) 

Fc=.02d0,Frk=l.d4,Frini=.02d0 

★ ■*■***★ ^r^^TUr******************************************-* ************ *★**★**★* 


Fc 

Frk 

Frini 


Tdepoc Epoch at initialization (days) modified to Julian 

date: 1950, Janl, 0hour=0.d0 

Tlnini Undeformed initially deployed tether length (km) 

Tlaini Deformed initially deployed tether length (km) 

Tlnini Undeformed tether length (km) initially stored in 

the subsatellite 

Tlnin2 Undeformed tether length (km) initially stored in 

the shuttle 

Dvnini Initial deploy/retract velocity from shuttle 

(km/sec) 

Dacini Initial deploy/ retract acceleration from shuttle 

(km/sec**2) 

Tdepoc=10400.d0,Tl ni ni=.ld-l,Tlaini=.ld-l 

Tlninl=0.d0,Tlnin2=19.990d0 

Dvnini=.2d-3,Daci ni=0.d0 


28 



Rshutt Cartesian position vector (km) of the shuttle 

referenced to the earth axis of 1950 
Vshutt Cartesian velocity vector (km/sec) of the shuttle 

referenced to the earth axis of 1950 
Tdirec Direction vector of the straight tether line 

referenced to the orbital reference frame with 
origin at shuttle, (automatically normalized) 
Omeini Angular velocity vector (1/sec) of the straight 

tether line referenced to the orbital reference 
frame 

Rshutt=6671 .d0,0.d0,0.d0 
Vshutt=0.d0,6.80d0,3.69d0 
Tdi rec=0.d0,0.d0, -1 ,d0 
0meini=0.d0,0.d0,0.d0 

Tedens Linear mass density of the tether (kg/km) 

Shmass Mass of the shuttle (kg) without tether 

Sumass Mass of the subsatellite (kg) without tether 

Tedens=.3d0,Shmass=900O0.d0,Sumass=18O.d0 

************************************************************************** 


Iteg 

Flag to set integration routine 

1 

Variable step size routine 

2 

4th order Runge-Kutta 

Iteg=l 


Nstop 

Stopping condition by step control 

<=0 

No stopping due to step control 

> 0 

Stops after Nstop steps 

Tstop 

Stopping condition by time control 

<=0 

No stopping due to time control 

> 0 

Stops after Tstop seconds 

Nprint 

Printing condition by step control 

<=0 

No printing due to step control 

> 0 

Printing after Nprint steps 

Tprint 

Printing condition by time control 

<=0 

No printing due to time control 

> 0 

Printing after Tprint seconds 


Nstop=8000,Tstop=10000.d0 

Nprint=20,Tprint=-100.d0 
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3.0 Tension Deployment 


3.1 Theory 

To implement tension controlled deployment as a new option in the tethered 
satellites simulation, it is first necessary to develop the theory and determine 
the program modifications necessary. It is most desirable that none of the cur- 
rent options of the program be disturbed and that simulation run speed should 
not be adversely affected. The original version of the simulation program used 
the assumption that the undeformed length of tether deployed at any time is 
known and is constrained to follow a user-defined function of time. As a con- 
sequence of this, the tether deployment acceleration couples extensively into 
the calculation of the time derivatives of the other dynamic variables. Since, 
the deployment acceleration is constrained, the tension at the deployer is 
determined and beyond the control of the user. To gain user control of the 

deployer tension, it is necessary that deployed, undeformed tether length be a 
dynamic variable. Thus, we must determine another dynamic equation for tether 
length. Before we do this let us review the dynamical equations of the simula- 
tion. 


Figure 3.1.1 shows the basic coordinate systems used in the simulation. 
The equations of motion are solved in the inertial reference frame. The X axis 
of this reference frame is directed toward the first point of Aries (approxi- 
mately the position of the sun as seen from the earth at the time of the Vernal 
Equinox). The Z axis is aligned with the earth spin axis and, thus, points to 
the celestial north pole. The Y axis is defined to complete a right handed 
coordinate system. Two other reference frames are used for output by the simu- 
lation. Subsatellite positions are defined relative to the shuttle in a local 
vertical reference frame defined as shown. This frame is oriented with the Z 
axis aligned with the local vertical position away from the earth. The Y axis 
points along the orbital angular velocity vector (normal to the orbit plane). 
The X axis completes the right handed set and points generally along the velo- 
city vector. The output components for the tether nodal vectors are given with 
respect to a special local vertical reference defined for consistency with 
current tether analysis convention. 


Simulation Coordinate Systems 


♦ Inertial Reference 

Z (Geographic 



• Local Vertical 



Y (Orbit Normal) 
X 


Z (Local Vertical) 


• Special Local Vertical 



Figure 3.1.1. 


These coordinate systems are the only ones of any consequence used by the simu- 
lation. 

Figure 3.1.2 depicts the geometry of the shuttle-tether-subsatellite sys- 
tem. The shuttle and subsatellite are treated as point masses. The tether is 
treated as an elastic continuum of nonzero mass. The individual mass elements 
of the tether are specified by a length(il) measured from the subsatellite. This 
parameter measures nominal distances along the tether and the value at any mass 
point corresponds to the unstretched arc length measured from the end. The 
deformed arc length along the tether will differ from (A) by the amount of 
stretch. Other parameters of interest are L, the total deployed (undeformed) 
length of tether and Lt, the total length of tether available for deployment. 
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Partial Differentia! Equation (PDE) 
Tethered Satellites Dynamics 
Model 


Shuttle 



L — Deployed tether length 

Lj — Total tether length 

— Mass point on tether 
measured from satellite 


Satellite 


Figure 3.1.2. 


32 



Figure 3.1.3 shows a typical element of tether at some position (1) along 
the tether. The position vector defining the position of (A) is r(l). The 
forces acting on this element are due to three sources. The first two are 
material forces due to tether elasticity: tension and shear. The third force 

is due to the external environmental forces acting on the element. These 
include earth gravity modelled by spherical and nonspherical earth terms, sun 
and moon gravity, aerodynamic drag and solar radiation pressure. The form of 
the tension and shear forces is shown in the figure. The other external forces 
are based on well-known models and are not discussed. They are explained in the 
original documentation of this simul ation.[l] 


N({+d!) 


Q(f+df) 



f(f) - Body force per unit length 
,Q(P) - Shear Force 


'■(p) " ^ —*-/ N(f) - (Normal Force) 

Tether Segment 


Tension: 

Shear: 


NW =/s[!r’(/)l -l]j^ ;l’ = 

' -1^ [Fi fej] - 


External Forces f(J) 


I ^ M 'p 


at 

oC 


m 


1. Spericoi Earth Gravity 

2. Non spherical Earth perturbations 

3. Sun & Moon Gravity 

4. Aerodynamic Dreg 

5. Soior Rodiation .pressure 


Figure 3.1.3. 

1. "Dynamics of a System of Two Satellites Connected by a Deployable and 

Extensible Tether of Finite Mass"; P. Kohler, W. Maag, and R. Wehrli, 1978. 
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Approximate solutions can be obtained by using appropriate numerical proce- 
dures. The partial differential equations appearing in Figure 3.1.4 are written 
in terms of two independent variables, time and undeformed length. The numeri- 
cal integration process determines solutions at discrete future times. The spa- 
tial variation must also be discretized to obtain the solution numerically. One 
approach to discretization would be to track the motion of mass points located 
along the tether at values that are multiples of some fixed distance. The dif- 
ficulty with this approach is that the number of such points (nodes) required 
varies with the length of tether deployed. This is inconvenient since the 
number of equations being integrated would change with the amount of tether 
deployed. To overcome this, we define a normalized length (Ci) = VL and trans- 
form the equations of motion to this variable. The new equations of motion are 
obtained by substitution of the results of the change of variables shown in 
Figure 3.1.5. If we now select nodal points to be located at the two ends and 
at points separated by 1/(N-1) for N nodes, we can write the equations of motion 
of the tether in terms of a constant number of nodes. This makes the job of 
numerical solution considerably simpler. On the other hand, many more terms are 
introduced into, the equations of motion to account for the variation of the 
length L. This accounts for the fact that a variable amount of mass is asso- 
ciated with a tether node depending on the length L. Defining new variables T 
and V allows us to express the equations of motion in first order form. The 
partial differential equations are converted to difference equations according 
to the procedures shown in Figure 3.1.6. This process results in the equations 
summarized in Figure 3.1.7 with auxilliary variables as defined on the bottom of 
Figure 3.1.7 and Figure 3.1.8. 

The set of equations summarized in Figures 3.1.7 and 3.1.8 represents the 
original form of the simulation as acquired by Control Dynamics under license 
from ACM/SAI. No material damping was included in this model. Test results 
suggest that energy dissipation in tether materials is significant. Thus, to 
model the material damping that is shown to be present, we added viscous damping 
to the tether tension model. This model is shown in Figure 3.1.9. 

A more difficult modification to the simulation was the addition of the 
capability to simulate tension controlled deployment. The simulation as origi- 
nally formulated treated the deployed length L as a known function of time 
specified by the user. From this known function and the dynamic equations for 
the system, the tension within each segment of the tether is determined. To 
make the tension controllable it is necessary to make the tether length a 
dynamic variable. This can be accomplished in several ways. The way chosen 
minimizes the required changes to the simulation and maintains all existing 
capabilities. The chosen technique is shown in Figure 3.1.10. 


The dynamic equations for the typical tether element are derived in Figure 
3.1.4 based on the vector diagram shown in Figure 3.1.3. The motion of the end- 
masses is defined by the boundary conditions. These boundary conditions are 
essentially just the dynamic equations for point masses. Body 2 (the shuttle) 
includes an extra term which looks like a thrusting term due to the flow effect 
of the deploying tether and results from the same sort of considerations used to 
derive the dynamical equations of a rocket and is referred to as the rocket term 
for this reason. If bending stiffness of the tether is nonzero, additional 
boundary conditions must be satisfied. These are shown at the bottom of Figure 
3.1.4. 


Dynamics Equations of MotioQ 
Boundary Conditions 

= N(f+dO - N(?) + O(i+d0 - Q(0 + f df 

= + Q(0] +1 df 

/^l(0 = + Q(J)j + f : Interior tether point 

Boundary Conditions: 
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=Ii+ [1(0 + Q(0J 

2 i' 2 = I 2 - Ki) -Q(f)] j 


L 


for«>< Y 0 


(r' 

(l' 




Jf = 0 ~ ° 
i= L = ° 


End body dynamics 


Bending moment must 
vanish at tether ends 


Figure 3.I.4. 
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DISCRETIZATION 

Figure 3.1.5 

Define normalized length 

Change to tether parameter from 1. 


r 0,t) = ? (^,t) 


r' 
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H const 



t const 



f = f + 2v'i + ?'C+ I'i 
Let T = r’ 

w _ jj Figure 3.1.6 


To discretize the spatial equations assume n uniformly 

distributed nodes including boundary nodes. 
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TENSION CONTROLLED DEPLOYMENT 


• Original Sim. Assumed L, L, L Constrained 


y - fo (y) + L f , (y) 
L == L(t) 

C = ^(t) 

C = L(t) 


Modified Sim. 


y = f. (y) + L f. (y) 

il' = C(y) + lL) y 


Fully Extensible, Homogeneous 


y = T(y) + L f,(y) 


F - N 


n- 1/2,0 


N 


n- 1/2,1 


Inextensible 


Figure 3.1.10 
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The details of the development of modifications to the equations of motion 
of tethered satellites simulation have been discussed in the previous para- 
graphs. These modifications have been implemented and verified. 


3.2 Implementation into Simulation 

The following additions and changes were made to the program for the 
constant tension deployment law. 

The vectors Y, Z, and F were expanded to allow for the integration of tether 
deployment acceleration and velocity and the friction force derivative. This 
change affects common UMESH and was changed in every subroutine in which these 
terms appear. This increased the dimensions of these variables to 123. 


COMMON/DEPLAW/ALPHA,IDEPL,FC,FRK,FR (added FC,FRK,FR) 

This common appears in the Main, Subroutine Deplex, Subroutine Rbegin, 
Subroutine Rinitl, Subroutine Rwrite, Subroutine Tendpl, Subroutine 
Uatox, Subroutine Ubegin, Subroutine Uder, Subroutine Ugen, Subroutine 
Uinitl, Subroutine Uprint, and Subroutine Uxtoa. 

Subroutine Deplex 

GO T0(10,20,30,40,50)IDEPL 
50 CONTINUE 

Subroutine Kinmax(FO,Fl) 

IMPLICIT REAL*8(A-H,0-Z) 

COMMON/SHUTCO/RSHUTT ( 3 ) , VSHUTT( 3 ) ,RSHU2 , RSHU 
C0MM0N/RELC0/TAU(19,3),RELVEL(19,3) 

COMMON/TDEPPA/TLN , DVN , DAC , XMASS2 
C0MM0N/TINDPA/DUMMY1(4),XINCRE,DUMMY2(2),N0DES 
COMMON/TESPEC/TEDENS, BETA, GAMMA, LONELA 
DIMENSION F(60),F0(60),F1(60) 

C 

NMINl = NODES-1 

NMIN2 = NODES-2 

N3=3*N0DES 

X2=DVN/TLN 

X3=-XINCRE*X2*X2 

DAC=0.D0 

1 XM=0.5D0*XINCRE*DAC/TLN 

X1=XM+X3 
XP=X1+X3 

C CONTRIBUTION TO THE DERIVATIVES OF RELVEL(1,I) 

DO 20 1=1,3 

F(I)=XM*TAU(1,I) + XP*TAU(2,I) + X2*(RELVEL(1,I)+RELVEL(2,I) ) 
20 CONTINUE 
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C CONTRIBUTION TO THE DERIVATIVES OF RELVEL(2, I ) , . . . ,RELVEL(NODES-2, I ) 

DO 30 K=2,NMIN2 
FL0AT1=DBLE(1-K) 

FL0AT2=DBLE(K) 

XX1=FLOAT1*(X1+FLOAT1*X3) 

XX2=XM+2 .D0*FL0AT1*FL0AT2*X3 
XX3=FL0AT2*(X1+FL0AT2*X3) 

DO 40 1=1,3 
J=3*(K-1)+1 

F { J )=XX1*TAU(K-1 , I )+XX2*TAU (K, I )+XX3*TAU(K+l , I ) 

$ + X2^(FL0AT1*RELVEL(K-1,I)+RELVEL{K,I)+FL0AT2*RELVEL(K+1,I)) 

40 CONTINUE 

30 CONTINUE 

C CONTRIBUTION TO THE DERIVATIVES OF RELVEL(N0DES-1,I ) 

FL0AT1=-DBLE(FL0AT(NMIN2)) 

XX1=FL0AT1*(X1+FL0AT1*X3) 

XX2=FL0AT1*(X1-FL0AT1*X3) 

DO 50 1=1,3 
J=3*NMIN2+I 

F(J)=XX1*TAU(NMIN2,I)+XX2*TAU(NMIN1,I)+FLA0T1*X2* 

$ (RELVEL(NMIN2,I)+RELVEL(NMIN1,D) 

50 CONTINUE 

C ROCKET TERM. CONTRIBUTES TO THE DERIVATIVES OF RELVEL(N0DES-1, 1 ) AND 

C VSHUTT(I) 

XX1=0.5D0*TEDENS*DVN*X2/XMASS2 
DO 60 1=1,3 

XX2 = XX2*(3.D0*TAU(NMIN1,I)-TAU(NMIN2,I)) 

J1 = 3*NMIN2+I 
J2 = Jl+3 
F(J1) = F(J1)+XX2 
F(J2) = XX2 
60 CONTINUE 

IF(DAC .LT. .5D0) THEN 
DO 70 I=1,3*N0DES 
FO(I)=F(I) 

70 CONTINUE 

DAC=1.0D0 
GOTO 1 
ENDIF 

DO 80 I=1,3*N0DES 
F1(I) = F1(I) - FO(I) 

80 CONTINUE 

RETURN 
END 
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Subroutine Rbegin 

READ(19,*,END=1) INSTEP, TDEPOC, TLNINI ,DVNINI .DACDUMJLNINl, 
$ TLNIN2,RSHUTT(1) ,RSHUTT(2) ,RSHUTT(3) ,VSHUTT(1) ,VSHUTT(2) , 

$ VSHUTT(3),FRINI 

G0T0(3,4,6,5,7),IDEPL 

3 CONTINUE 

4 DACINI=0.D0 
GOTO 6 

5 CONTINUE 
DACINI=ALPHA*DVNINI 
GO TO 6 

7 TLN=TLNINI 

DVN=DVNINI 
DAC=DACDUM 
FR=FRINI 

6 CONTINUE 

WRITE(16,*) 'FC: ' ,FC, *FRK: ' ,FRK 

Subroutine Riniti 

READ(5,*)FC,FRK,FRINI 

PRINT *, 'FC,FRK,FRINI; ' ,FC,FRK,FRINI 

FC = FC/IOOO.DO 

FRINI = FRINI/IOOO.DO 

FRK = FRK/IOOO.DO 


Subroutine Rwrite 

WRITE(19,*) NSTEP, TDNOW, TIN, DVN, DAC, TLNl, TLN2, 

$ RSHUTT(l), RSHUTT(2), RSHUTT(3), VSHUTT(l), VSHUTT(2), 

$ VSHUTT(3), FR 

Subroutine Tendpl (T,XL0DD,XL10D) 

THIS SUBROUTINE COMPUTES TIME DEPENDENT PARTS FOR A 
TENSION DEPLOYMENT LAW. THE PARAMETERS TO BE COMPUTED ARE 
DVN: CONSTRAINED DEPLOYMENT RATE BASED ON TENSION 

XLODD: TETHER ACCELERATION EXCLUDING TERMS IN Z 

XLlDD(l-3): COEFFICIENTS OF RPDD AT NODE N 

IMPLICIT REAL*8(A-H,0-Z) 

DIMENSION RP(3),RPP(3),RPD(3),RPPD(3),XL1DD(3) 
COMMON/TINDPA/TINI, TLNINI, DVNINI,DACINI,XINCRE,XMASSO,XMASSl, 
$ NODES 

COMMON/TDEPPA/TLN , DVN , DAC , XMASS2 
COMMON/DEPLAW/ALPHA,IDEPL,FC,FRK,FR 
C0MM0N/RETRAP/TLNIN1,TLNIN2,TLN1,TLN2 
COMMON/EMASS/ SUMASS.SHMASS 
C0MM0N/NF0RCE/XNF0RC(19) 

COMMON/INTERM/TAUQUA(19) ,TAUMA6(19) ,BLABLA(37) 
COMMON/UMESH/TT,DTT,Y(123),F{123) 
C0MM0N/SHUTC0/RSHUTT(3),VSHUTT(3),RSHU2,RSHU 
C0MM0N/RELC0/TAU(19,3),RELVEL(19,3) 
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c 

NMINl = NODES - 1 

NMIN2 = NMINl - 1 

NMIN3 = NMIN2 - 1 

XNMINl = DBLE(NMINl) 

C COMPUTE AUXJLLIARY VARIABLES 

XI = l.DO - .5D0*XINCRE 
C COMPUTE MORE AUXILLIARY VARIABLES 

DO 10 1=1,3 
RP(I) = TAU(NMINl.I) 

RPD(I) = RELVEL(NMIN1,I)/XINCRE 

RPP(I) = (1.5D0*TAU(NMIN1,I) - 2.0D0*TAU(NMIN2,I) + 

$ .50D0*TAU(NMIN3,I))/XINCRE 

10 RPPD(I) = (1.5D0*RELVEL(NMIN2,I) - 2.0D0*RELVEL(NMiN2,I) + 

$ .5D0*RELVEL(NMIN3( I ) )/XINCRE/XINCRE 

RPMSQ = RP(1)**2 + RP(2)**2 + RP(3)**2 
RPM = DSQRT(RPMSQ) 

C COMPUTE LOOT NUMERATOR TERMS 

XNl = TLN*RPMSQ 
XN2 = RPM*TLN*TLN 

XN3 = TLN*(RP(1)*RPD(1) + RP(2)*RPD(2) + RP(3)*RPD(3) ) 

XN4 = FR*XN2 

C COMPUTE LOOT DENOMINATOR TERMS 

XDl = XI*(RP(1)*RPP(1) + RP(2)*RPP(2) + RP(3)*RPP(3) ) 

XD2 = RPMSQ 

XNUM = (BETA*(XN1-XN2) + GAMMA*XN3 - XN4) 

XDENOM = GAMMA* (XDl + XD2) 

C COMPUTE LOOT (DVN) FROM FR 

DVN = XNUM/XDENOM 
IF (DABS(FR) .LT. FC) THEN 
FRDOT = FRK*DVN*TAUMAG(NMIN1)/TLN 
ELSE 

FRDOT = O.DO 
ENDIF 

C DERIVATIVES OF NUMERATOR AND DENOMINATOR TERMS 

XNID = XN1*DVN/TLN + 2.D0*XN3 
XM2D = 2.D0*DVN*XN2/TLN + TLN**3*XN3/XN2 
XN3D = DVN*XN3/TLN 
XDID = O.DO 
XD2D = 2.D0*XN3/TLN 
DO 20 1=1,3 

XN3D = XN3D + TLN*RPD(I )*RPD( I ) 

XDID = XDID + XI*RPP(I)*RPD(I) + XI*RP( I )*RPPD(I ) 

XLIDD(I) = TLN*RP(I)/(XD1 + XD2) 

20 CONTINUE 

XN4D = FR*XN2D + FRD0T*XN2 

XNUMD = BETA*(XN1D-XN2D) + GAMMA*XN3D - XN4D 

XDEND = GAMMA* (XDID + XD2D) 

XLODD = (XNUMD-XNUM*XDEND/XDENOM)/XDENOM 

RETURN 

END 
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Subroutine Uatox 

IF(I0EPL .EQ. 5) THEN 
FR = Y(6*N0DES+1) 

TIN = Y(6*N0DES+2) 

IF(L0NELA .EQ. 1) DVN = Y(6*N0DES+3) 

IF(DABS(FR) .GT. FC) FR = SI6N(FC,FR) 

ENDIF 

Subroutine Ubegin 

G0T0(3,4,6,5,7)IDEPL 

3 DVNINI=O.DO 

4 DACINI=0.D0 
GOTO 6 

5 DVNINI=ALPHA*TLNINI 
DACINI=ALPHA*DVNINI 
GO TO 6 

7 TLN = TLNINI 

DVN = DVNINI 
DAC = DACINI 

6 CONTINUE 

WRITE(16,*)'FC: ‘,FC,'FRK: ',FRK 

Subroutine Uder 

COMMON/DEPLAW/ALPHA.IDEPL,FC,FRK,FR 
COMMON/TDEPPA/TLN , DVN , DAC , XMASS2 
C0MM0N/RETRAP/TLNIN1,TLNIN2,TLN1,TLN2 
C0MM0N/INTERM/TAUQUA(19),TAUMAG(19),TIIP1(18),RELVQU(19) 
COMMON/NFORCE/XNFORC (19) 

DIMENSION XN0(19),XN1(19) 

NMIN2 = NODES - 2 
NMIN3 = NODES - 3 
CALL HELP 

IF(IDEPL .EQ. 5 .AND. LONELA .GT. 1) THEN 
CALL TENDPL(T,XL0DD,XL1DD) 

ENDIF 

C COMPUTE NORMAL FORCE FOR LDD (DAC) = 0 

DAC = O.DO 
DO 260 0=1, N3 
JP=J+N3 

Z(JP)=Z(J)+F0(J) 

FP(0)=Z(P) 

260 CONTINUE 

CALL NORMAL (Z,F) 

DO 265 0=1, N3 
0P=0+N3 

Z(0P)=Z(0P)+F(0) 

F0(0)=Z(0P) 

265 CONTINUE 

DO 270 I=1,NMIN1 
XNO(I)=XNFORC(I) 

270 CONTINUE 
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c 

280 

285 

290 

1005 

295 

300 

C 

C 

40 

C 

310 


COMPUTE NORMAL FORCE FOR LOO (OAC) = 1 

DAC=1.D0 

DO 280 J=1,N3 

JP=J+N3 

Z(JP)=FP(J)+F1(J) 

FP(J)=Z(JP) 

CONTINUE 
CALL NORMAL (Z,F) 

DO 285 J=1,N3 
JP=J+N3 

2(JP)=Z(JP)+F(J) 

F1(J)=Z(JP) 

CONTINUE 

DO 290 I=1,NMIN1 
XN1(I)=XNF0RC(I) 

CONTINUE 

IF(LONELA .EQ. 1) THEN 

DAC=(FR/TAUMAG(NMIN1)-XN0(NMIN1))/(XN1(NMIN1)-XN0(NMIN1)) 

ELSE 

DAC=XL0DD 
XDEN=1.D0 
DO 1005 1=1,3 

XDEN=XDEN-XL1DD( I )*(F1 (N3-6+I )-F0(N3-6+I ) )/XINCRE 

DAC=DAC+XLlDD(I)*F0(N3-6+I)/XINCRE 

CONTINUE 

DAC=DAC/XDEN 

ENDIF 

DO 295 J=1,N3 
JP=J+N3 

Z(JP)=FO(J)+DAC*(F1(0)-FO(J)) 

CONTINUE 

DO 300 I=1,NMIN1 

XNFORC( I )=XNO( I )+DAC*(XNl ( I )-XNO( I ) ) 

CONTINUE 

INTEGRATION VARIABLES FOR TENSION DEPLOYMENT 
Z(6*N0DES+1) = FRK*DVN*TAUMAG(NMIN1)/TLN 
Z(6*N0DES+2) = DVN 

IF(LONELA .EQ. 1) Z(6*N0DES+3) = DAC 
ELSE 

COMPUTE KINEMATICAL TERMS 
CALL KINEMA(F) 

DO 40 J=1,N3 
JP=J+N3 

Z(JP)=Z(JP)+F(J) 

CONTINUE 

THE NORMAL FORCE 
CALL NORMAL (Z,F) 

DO 310 J=1,N3 
JP=J+N3 

Z(JP)=Z(JP)+F(J) 

CONTINUE 

ENDIF 
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Subroutine Ugen 

IF(I0EPL .EQ. 5) THEN 
NDEQ = 6*N0DES + 2 

IF(L0NELA .EQ. 1) NDEQ = 6*N0DES + 3 
ENDIF 


Subroutine Uiniti 

READ(5,*)FC,FRK,FRINI 

PRINT *,'FC,FRK,FRINI: ',FC,FRK,FRINI 

FC = FC/IOOO.DO 

FRINI = FRINI/IOOO.DO 

FRK = FRK/IOOO.DO 

IF(IDEPL .EQ. 5) FR = FRINI 

Subroutine Uxtoa 

IF(IOEPL .EQ. 5) THEN 
Y(6*N0DES+1) = FR 
Y(6*N0DES+2) = TLN 
Y(6*N0DES+3) = DVN 
ENDIF 
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4.0 Damping at Deployment Point 

4.1 Theory 

While working on the tension deployment scheme, it became apparent that 
some type of damping would be necessary to reduce the effects of the high fre- 
quency modes travelling through the tether. The high frequencies led to extre- 
mely small time steps and long computer run times. The damping keeps the basic 
tether motions while reducing high frequency effects. 


The damping terms were incorporated in Subroutine Fulext for nodes 1 
through n-1. Subroutine Homoge for node n-1 (current coding then applies damp- 
ing to the remaining nodes). As no stretch is obtained in the inextensible 
tether, it was not necessary to add damping. 


To obtain the damping equation, it was first 
Rayleigh Dissipation Function, F. This term is similar 
energy equation with 3 replaced with Y and derivatives 
the equation 



necessary to generate a 
in form to the potential 
of used. This yields 


Take the partial with respect to £' to obtain the damping terms to incorporate 
in the force equations. 


8F 

3r‘ 


Y 

L 


CL 

[ _T ' ._T + _T.2 

L 



T 



2 


In these equations: 


L -*■ deployed tether length 
L tether deployment velocity 
Y damping coefficient 
C * floating tether point 
2 position of floating tether point. 
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4,2 Implementation into Simulation 

The following additions and changes were made to the program for the 
tether damping. 

COMMON/TESPEC/TEDENS, BETA, GAMMA, LONELA (GAMMA added in) 

This change occured in the MAIN, Subroutine Deplex, Subroutine 
Fulext, Subroutine Head, Subroutine Help, Subroutine Homoge, Subroutine 
Kinema, Subroutine Kinmax, Subroutine Ngrav, Subroutine Normal, Subroutine 
Rbegin, Subroutine Rinitl, Subroutine Tendpl , Subroutine Uatox, Subroutine 
Ubegin, Subroutine Uder, Subroutine Ugen, Subroutine Uinitl, and Subroutine 
Upri nt. 


Subroutine Tendpl has previously been entered in the tension 
deployment section of the report, the damping terms were included. 


Subroutine Fulext 

C0MM0N/RELC0/TAU(19,3),RELVEL(19,3) 

NMIN2 = NODES-2 
NMIN3 = NODES-3 
DO 10 K=1,NMIN1 
X=BETA*(X1-1.D0/TAUMAG(K) ) 

XD = RELVEL(K,1)*TAU(K,1) + RELVEL(K,2)*TAU(K,2) + RELVEL(K,3)* 
$ TAU(K,3) 

IF(K .LT. NMINl) THEN 
IF(K .EQ. 1) THEN 

TAUSQP = .25D0*(-3.D0*TAUQUA(K) + 4. DO* 

$ TAUQUA ( K+1 ) -TAUQUA ( K+2 ) ) /X I NCRE 

ELSE 

TAUSQP = .25D0*(TAUQUA(K+1) - TAUQUA(K-l)) 

$ /XINCRE 

ENDIF 
ELSE 

TAUSQP = .25D0*(3.D0*TAUQUA(NMIN1) - 4. DO* 

$ TAUQUA(NMIN2) + TAUQUA ( NMI N3 )) /XI NCRE 

ENDIF 

XD = XD/XINCRE - DVN/TLN*(TAUQUA(K) + DBLE(FL0AT(K) ) - .5D0)* 

$ TAUSQP/DBLE(FL0AT(NMIN1))) 

XD =• GAMMA*XD/TLN/TAUQUA(K) 

X = X + XD 

IF(ISLACK .EQ. 2 .AND. X .LT. O.DO) X=0.D0 
XNFORC(K) = X 
10 CONTINUE 

Subroutine Homoge 

C0MM0N/RELC0/TAU(19,3),RELVEL(19,3) 

NMIN3 = NODES-3 

CC ADD DAMPING TO SECTION NEXT TO BODY 2 

C 

XD = RELVEL(NMIN1,1)*TAU(NMIN1,1) + RELVEL(NMIN1,2)* 

$ TAU(NMIN1,2) + RELVEL(NMIN1,3)*TAU(NMIN1,3) 

TAUSQP = .2500*(3.D0*TAUQUA(NMIN1) - 4.D0*TAUQUA(NMIN2). + 
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$ TAUQUA(NMIN3))/XINCRE 

XD = XD/XINCRE - DVN/TLN*(TAUQUA(NMIN1) + {l.DO - 
$ .5D0/DBLE(NMIN1))*TAUSQP) 

XD = GAMMA*XD/TLN/TAUQUA(NMIN1) 

Subroutine Rbegin 

READ ( 1 9, * ) LONNEW , TEDENS , BST I FF , BETA , GAMMA, NODES , SHMASS , SUMASS 
WRITE(16,*) 'LONELA: ' ,LONELA, 'BETA: ' .BETA, 'GAMMA* , GAMMA, 

$ * ISLACKt * ISLACK 

WRITE(19,*)L0NELA,TEDENS,BSTIFF, BETA, GAMMA, NODES, SHMASS, SUMASS 

Subroutine Rinitl 

READ(5,*)L0NELA, BETA, GAMMA, ISLACK 
PRINT *, LONELA, BETA, GAMMA, ISLACK 
GAMMA = GAMMA/ 1000. DO 

READ(19,*)L0NNEW,TEDENS,8STIFF, BETA, GAMMA, NODES, SHMASS, SUMASS 
Subroutine Ubegin 

WRITE(16,*) 'LONELA: ' .LONELA, 'BETA: ' .BETA, 'GAMMA: ' .GAMMA, 

$ 'ISLACK: '.ISLACK 

WRITE(19,*)L0NELA,TEDENS,BSTIFF, BETA, GAMMA, NODES, SHMASS, SUMASS 

Subroutine Uinitl 

READ(5,*)L0NELA, BETA, GAMMA, ISLACK 

PRINT *, 'LONELA, BETA, GAMMA, ISLACK: '.LONELA, BETA, GAMMA, ISLACK 
GAMMA = GAMMA/ 1000. DO 



5.0 Tektronics Graphics 


With the simulation running on the'VAX at MSFC, it became necessary to 
implement some type of graphics compatible with the Tektronics terminals. A 
basic scheme was first implemented so that the most fundamental types of graphs 
could be viewed. Throughout the project this initial code was revised and 
updated to its present form. (See Appendix-D) The graphics now include such 
features as a simple menu which returns to the screen after each plot, the 
ability to set limits on each axis for a plot, the ability to enter labels for 
the plots, and the ability to plot 'walking-plots'. The 'walking plots' are an 
attempt to emulate the plots yielded by the Smithsonian simulation. They depict 
the inplane/out-of-plane motion of the tether as it deploys. The ‘walking- 
plots' include additional features which allow the user to look at every nth 
point or scale the data as desired. This last option helps the user to better 
see the deflections in the tether as it deploys. 


5.1 Examples 

The following plots and graphs are examples of the graphics capabilities 
currently in use at MSEC on their VAX in coordination with the tether simula- 
tion. When running the plotting program, the first thing the user sees is a 
very basic menu (Figure 5.2.1) which lists the variables available and asks what 
type of plot the user desires. The user can specify either a normal x vs y type 
plot with up to four y variables, the 'walking-plots', a plot of tension at all 
the tether nodes (disregarding the shuttle as a node), or no plot at all. 

If the user requests the x vs. y type plot, the next question asks how 
many plots per graph, up to four, are desired. If a number not one through four 
is input, the program returns to the basic menu. The user inputs the plot per 
graph number. The program then asks for the x-axis number associated with the x 
variable desired, and then the corresponding y-axis numbers. Eight character 
labels for the x and y variables are requested next (Figure 5.2.2), and finally 
the program gives the current plot limits and asks if any change is desired. If 
so, then the user inputs the new plot limits. For one plot per graph, new 
limits are not requested. The plot then appears on the screen, is sent to the 
thermal printer, disappears, and the menu is returned. Any of the variables 
listed in the basic menu may be plotted in this manner. (Figure 5.2.3, tether 
length in km vs. time in sec.; Figure 5.2.4 in-plane and out-of-plane angles in 
degrees vs. time in sec.; Figure 5.2.5 tether deployment velocity in km/sec vs. 
time in sec.; Figure 5.2.6 tether deployment acceleration in km/sec**2 vs. time 
in sec.; Figure 5.2.7 the distance between the shuttle and the subsattelite in 
km vs. time in sec.; and Figure 5.2.8 the friction force applied at the 
deployment point in km vs. time in sec.) 

For the 'walking-plots' the user needs to know the number of nodes used 
in the simulation run as this is the first information requested when this type 
of plot is chosen. The next question asks if the user wants to view a plot of 
every time point, if not then a number n is requested to look at every nth data 
set. The first time through the 'walking-plot' section, the code calculates the 
maximum in-plane and out-of-plane deflections in the tether and where these 
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deflections occur. The user can then enter scale factors based on these deflec- 
tions. The user then specifies an in-plane or an out-of-plane graph. As 
before, the data is plotted to the screen, "sent to the printer, disappears, and 
the main menu returns. Figure 5.2,9 gives the questions asked when a walking 
plot is desired. Figure 5.2.10 is the plot associated with the above input, 
this is the basic plot with the scale factors at magnitude one and every time 
point plotted. Figure 5.2.11 shows this plot with every third time point 
plotted. Figure 5.2.12 has all the time points plotted, but the scale factor is 
.075 versus the previous value of 1.00. 

When a tension plot is requested, the program first asks for the base 
number of variables; i.e. the number of variables above the tension variables in 
the menu. In the example base menu this number would be 15. The code then 
plots the remaining tension variables on one graph. When finished, this process 
also returns to the main menu. Figure 5.2.13 shows the tension for nodes 1 

through the node next to the shuttle. These plots were obtained using the 
sample input file previously given. The markedly different behavior exhibited 
at the end of the run is accounted for in that the tether became fully deployed 
at approximately 9100 seconds. At this point the tether switched into the 

constant length mode for the duration of the run. 
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Figure 5.2.13 


6.0 Conclusions 


This report has presented the changes and modifications implemented to the 
tether simulation to form Version 3.0 as currently in existence at MSFC. 

As more changes are made to the simulation at MSFC, they will be documented 
and submitted as attachments to this report. Other updates to be included at a 
later time include the SEDS analysis cases, the AMIGA graphics task, the 
deployer mechanism simulation, and the spinning tether analysis. 
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APPENDIX A. SIMULATION CORRECTIONS 


Subroutine N6RAV: 

Original code 

DO 100 J=l,3 

TAUIMl(J) = TAU(J) 

TAU(J) = SNGL(FTAU(I,J)) 

TAUI(J) = (TAU(J) + TAUIMl(J))/2. 
LENI = LENT + TAUI ( J )*TAUI (J ) 

POS(J) = SNGL(FRNOD(I,J))*1000. 

100 CONTINUE 

Corrected code 

DO 100 JJ=1,3 

TAUIMl(JJ) = TAU(JJ) 

TAU(OJ) = SNGL(FTAU(I,JJ)) 

TAUI(JJ) = (TAU(JJ) + TAUIMl(JJ))/2. 
LENI = LENI + TAUI(JJ)*TAUI(JJ) ‘ 
POS(JJ) = SNGL(FRN0D(I,JJ))*1000. 

100 CONTINUE 

Original code 

VIT(J) = SNGL(FVN0D(I,J))*1000. 
Corrected code 

DO 61 JJ=1,3 

VIT(JJ) = SNGL(FVN0D(I,J0))*1000. 

61 CONTINUE 

Original code 

VN(L) = V(L) - VTAU*TAU(L)/TAUTAU 
Corrected code 

VN(L) = V(L) - VTAU*(TAU(L)/TAUTAU) 
Original code 

DO 250 L=l,3 

250 VN(L) = VNMAG*VN(L) 

Corrected code 

removed from code 


The variable J is being 
used already in the 
subroutine. 


The variable J is being 
used already in the 
subroutine. 


To avoid math problems. 


Redundant step. 


Original code Incorrect term. 

DO 40 L=l,3 

LL = (I-l)*3 + L 

IF(IDR) F(LL) = F(LL) + DBLE(-.5*RHO*CN*AN*VNA*VN(L)/(B*1000. ) ) 
IF(IRA .AND. .NOT. lECL) F(LL) = F(LL) + DBLE(R(L)/(B*1000. ) ) 

40 CONTINUE 
Corrected code 


DO 40 L=l,3 

LL = (I-l)*3 + L 

IF(IDR) F(LL) = F(LL) + DBLE(-.5*RH0*CN*AN*VNMAG*VN(L)/(B*1000. ) ) 
IF(IRA .AND. .NOT. lECL) F(LL) = F(LL) + DBLE(R(L)/(B*1000.) ) 

40 CONTINUE 

Original code Mixed terms. 

AN = 2*LENI*SNGL(FTERAD) 

Corrected code 

AN = 2.*LENI*SNGL(FTERAD) 

Subroutines Uinitl and Rinitl 
Original code 

not originally there 


Incorrect units. 



Corrected code 

BETA = RETA/IOOO.DO 
Subroutines Uprint and Trackf 

Original code Incorrect dimensions 

COMMON/ I NTERM/TAUQUA (19), TAUMAG (18), BLABLA ( 37 ) 

Corrected code 

COMMON/ I NTERM/TAUQUA (19), TAUMAG (19), BLABLA ( 37 ) 

Subroutine Uthird 

Original code Incorrect dimensions 

DIMENSION P(l) 

Corrected code 

DIMENSION P(7) 



APPENDIX B. SIMULATION DELETIONS 


Main: 


COMMON/PRAZIS/PREC 

C0MM0N/TINDPA/TINI,TLNINI,DVNINI,DAC1NI,XINCRE,XMASS1, NODES 
COMMON/CPOEMS/XMU , WMU , WMU I N , XMUM, XMUS 
COMMON/CONTRO/ICONT 
COMMON/CGEOME/RE,FLAT,FLATSQ,CLIGHT 

COMMON/CSOLAR/SOLPR,ECL,CECL,SECL,XLSUNA,XLSUNB,ECLECC,ECLOM 

C0MM0N/C0NTRL/CNTRL(3) 

C0MM0N/CDYNAE/STD50,0MT50,0MQ50,0MR0T,STD50R,0MT50R,0MQ50R 
COMMON/UMESH/T ,DT , Y ( 992 ) , Z ( 124 ) 

PREC = O.ODO 
Subroutine Airden: 

C0MM0N/JP0SMS/FILL(36),SUNV(4),LUNVEC(4) 

Subroutine Bendin 

C0MM0N/TINDPA/DUMMY1(4),XINCRE,DUMMY2,N0DES 
Subroutine Centra 

COMMON/CPOEMS/XMU .DUMMY ( 4 ) 

C0MM0N/T1NDPA/DUMMY1(4),XINCRE,DUMMY2,N0DES 
TYPE *,'IN CENTRA: ' ,RNODEA(K)=‘ 

TYPE *,K,RN0DEA(K) 

Subroutine Coot: 

C0MM0N/CGE0ME/RE,FLAT,FLATSQ,CLI6HT 

C0MM0N/CDYNAE/STD50 , 0MT50 , 0MQ50 , OMROT , STD50R , 0MT50R , 0MQ50R 
COMMON/CPOEMS/XMU , WMU , WMU I N , XMUM , XMUS 

COMMON/CSOLAR/SOLPR,ECL,CECL,SECL,XLSUNA,XLSUNB,ECLECC,ECLOM 

RE = MEAN RADIUS OF EQUATOR; FLAT = FLATTENING COEFFICIENT 

REF: GODDARD EARTH MODEL 6 

RE = 6378.144D0 

RE = 6378.156D0 

FLAT = 1. DO/298. 257D0 

FLATSQ = FLAT* (2. DO-FLAT) 

CLIGHT = SPEED OF LIGHT. REF: SAO SPEC. REP. 353, 1973. 

CLIGHT = 299792. 5D0 

STD50, 0MT50, 0MQ50 GIVE THE ANGLE BETWEEN GREENWICH MERIDIAN 

AND MEAN EQUINOX (DEG) FROMD = MODIFIED JULIAN DATE (1950) 
ACCORDING TOW = STD50 + 0MT50*D + 0MQ50*D*D 

STD50R, 0MT50R, 0MQ50R ARE THE SAME ANGLE EXPRESSED IN RADIANS 
REF: THE ASTRONOMICAL EPHEMERIS 1976, P. 531. 

STD50 = 100.075542d0 
0MT50 = 360.985647335D0 
0MQ50 = .29D-12 
STD50R = STD50*RAD 
0MT50R = 0MT50*RAD 
0MQ50R = 0MQ50*RAD 

OMROT = MEAN VELOCITY OF ROTATION OF THE EARTH (RAD/SEC) 

OMROT = 0MT50R/864.D2 

XMU = CENTRAL EARTH POTENTIAL. REF: SMITHSONIAN STANDARD EARTH III 

AND GODDARD EARTH MODELS 5 AND 6 
XMU = 398601. 3D0 
WMU = DSQRT(XMU) 

WMUIN = l.DO/WMU 



XMUM, XMUS = MOON, SUN POTENTIAL. REF: JPL DEV. EPHEMERIOES NO. 19 
XMUM = XMU/81.301D0 
XMUS = 132715. OD6 

SOLPR=RADIATION PRESSURE AT MEAN EARTH DISTANCE FROM THE SUN 
IN KG*KM/(M*SEC)**2 REF: NASA TM-X 64627 
SOLPR = 4.51D-9 

ECL, CECL, SECL= INCLINATION OF THE ECLIPTIC AND ITS COS AND SIN. 
XLSUNA, XLSUNB, ECLECC, ECLOM GIVE WITH A PRECISION OF 0.01 DEG 
A VALUE FOR THE LONGITUDE OF THE SUN (LS{RAD)) IN THE ECLIPTI 
FROM: LM=XLSUNA+XLSUNB*DAY; LS=LM+ECLECC*DSIN(LM-ECLOM) 

REF: THE ASTR. EPHEM. 1976, EXTRAPOLATED TO MOD 10000 

ECL = RAD*23.442224D0 

ECL = 0.4090619414299053D0 

CECL = DCOS(ECL) 

SECL = DSIN(ECL) 

XLSUNA = RAD*280.08120D0 
XLSUNB = RAD*.9856473389D0 
ECLECC = 3.343724E-2 
ECLOM = RAD*282.55137D0 
Subroutine Deplex 

C0MM0N/TINDPA/TINI,TLNINI,DVNINI,DAVINI,XINCRE,XMASS1, NODES 
Subroutine Fulext 

C0MM0N/TINDPA/DUMMY1(6),N0DES 
Subroutine Ghalfi 

COMMON/TI NDPA/DUMMYl ( 4 ) , X I NCRE , DUMMY2, NODES 
Subroutine Ginteg 

COMMON/T I NDPA/DUMMY 1 ( 4 ) , X I NCRE , DUMMY2 , NODES 
Subroutine Head 

C0MM0N/TINDPA/TINI,TLNINI,DVNINI,DACINI,XINCRE,XMASS1, NODES 
Subroutine Help 

C0MM0N/TINDPA/DUMMY1(4),XINCRE,XMASS1,N00ES 
Subroutine Honioge 

C0MM0N/TINDPA/DUMMY1(4),XINCRE,DUMMY2, NODES 
Subroutine Inext 

C0MM0N/TINDPA/DUMMY1(4),XINCRE,DUMMY2,N0DES 
Subroutine Kinema 

C0MM0N/TINDPA/DUMMY1(4),XINCRE,DUMMY2,N0DES 
Subroutine Ngrav 

COMMON/PRAZIS/PREC 

C0MM0N/JP0SMS/FILL(36),FFS0(4),FFLU(4) 

PRE' = SNGL(PREC) 

DO 500 1=1,19 

500 TYPE *,'IN NGRAV, FRNOA( ' ,1, ' )=' ,FRNOA(I) 

Subroutine Normal 

COMMON/T I NDPA/DUMMY ( 4 ) , X I NCRE , XMASSl , NODES 
DIMENSION Z(124),F(60) 

Subroutine Pertur 

C0MM0M/TINDPA/DUMMY(7), NODES 
TYPE *,'INPERTUR: K,RNODEA(K)‘ 

TYPE *,K,RNODEA(K) 


Subroutine Radial 

COMMON/PRAZIS/PREC 

WR ITE ( 6 , 1 100 ) AERA , REFCO , CAL , TRANSFVSCON 
PRE = SNGL(PREC) 

IF(ABS(CAL) .LE. PRE) CAL = 0.0 
Subroutine Rbegin 

C0MM0N/TINDPA/TINI,TLNINI,DVNINI,DACINI,XINCRE,XMASS1, NODES 
C0MM0M/C0NTRL/CNTRL(3) 

REA0(5,*)CNTRL 

$ VSHUTT ( 3 ) , CNTRL ( 1 ) , CNTRL ( 2 ) , CNTRL ( 3 ) 

Subroutine Riniti 

C0MM0N/TINDPA/TINI,TLNINI,DVNINI,DACINI,XINCRE,XMASS1, NODES 

COMMON/CPOEMS/XMU,WMU,WMUIN,XMUM,XMUS 

COMMON/CONTRO/ICONT 

IC0NT=1 

Z0NAL(2) = O.DO 
Subroutine RK78 

SUBROUTINE RK78( IR,T,DT,X,DUM,F1,F2.F3,F4,F5,F6,F7,N, 

$ TOL.DER) 

DIMENSION X(124),DUM(124),F1(124),F2(124),F39124),F4(124), 

$ F5(124),F6(124),F7(124),T0L(124),CH(13),ALF(10) 

PRINT *,'X7 (MUST BE <= 1 TO PASS DT),DT' ,X7,DT 
PRINT *,'DT ACCEPTED; DT,T:',X9,T 
Subroutine ROT 

COMMON/PRAZIS/PREC 
PRE = SNGL(PREC) 

Subroutine Rwrite 

COMMON/CONTRL/CNTRL(3) 

COMMON/TINDPA/DUMMYl ( 4 ) , XINCRE , XMASSl , NODES 
$ VSHUTT(3),CNTRL(1),CNTRL(2),CNTRL(3) 

Subroutine Thrust 

COMMON/TINDPA/DUMMY ( 4 ) , X I NCRE , XMASSl , NODES 
Subroutine Trackf 

COMMON/UMESH/T , DT , Y ( 992 ) , Z ( 1 24 ) 

COMMON/TINDPA/TINI,TLNINI,DVNINI,DACINI, XINCRE, XMASSl, NODES 
Subroutine Uatox 

COMMON/TINDPA/DUMMYl (4), XINCRE, DUMMY2, NODES 
DIMENSION Y(124) 

TYPE *, 'IN UATOX: FOR I = ',I 

TYPE *,'RN0DE(N-1,I)=',RN0DE(NMIN1,I), 

$ 'VN0DE(N-1,I)=',VN0DE(NMIN1,I) 

TYPE *,'RNODE(',L,',I)=‘,RNODE(L,I), 

$ ‘VNODE(',L, VN0DE(L, I) 

TYPE *,'IN UATOX: K,RNODEA=‘ 

TYPE *,K,RNODEA(K) 

Subroutine Ubegin 

COMMON/TINDPA/TINI,TLNINI,DVNINI,DACINI, XINCRE, XMASSl, NODES 
C0MM0N/C0NTR0L/CNTRL(3) 

READ(5,*)CNTRL 
WRITE(16,*) 'CNTRL: ', CNTRL 



•Subroutine Uder 

C0MM0N/TINDPA/TINI,TLNINI,DVNINI,DACINI,XINCRE,XMASS1, NODES 

COMMON/PRAXIS/PREC 

DIMENSION Y(124),Z(124),F(60) 

DO 15 K=1,NMIN1 
DO 15 1=1,3 
X=DABS(TAU(K,I)) 

IF(X.LT.PREC)TAU(K,I)=O.DO 

X=DABS(RELVEL(K,I)) 

IF(X.LT.PREC)RELVEL(K,I)=O.DO 
15 CONTINUE 
Subroutine Ugen 

C0MM0N/UMESH/T,DT,Y(1116) 

COMMON/TINDPA/TINI,TLNINI,DVNINI,DACINI,XINCRE,XMASSl, NODES 
Subroutine Uinitl 

C0MM0N/TINDPA/DUMMY1(4),DUMMY2(2),N0DES 
COMMON/CPOEMS/XMU , WMU , WMU I N , XMUM , XMUS 
COMMON/CONTRO/ICQNT 
ICONT=l 

ZONAL(2) = O.DO 
Subroutine Uint 

DIMENSION TOL(124) 

C0MM0N/UMESH/T,DT,A(124),F1(124),F2(124),F3(124),F4(124), 

$ F5(124),F6(124),F7(124),DUM(124) 

CALL RK78(IR,T,DT,A,DUM,F1,F2,F3,F4,F5,F6,F7,NDEQ,T0L,UDER) 
Subroutine Unosph 

COMMON/CD YNAE /STD50 , 0MT50 , 0MQ50 , OMROT , STD50R , 0MT50R , 0MQ50R 
Subroutine Uprint 

C0MM0N/UMESH/T,DT,Y(992),Z(124) 

COMMON/TINDPA ( DUMMYl (4 ) , X I NCRE , XMASSl , NODES 
TYPE *,'NSTEP: ',NSTEP,’TIME',T 
Subroutine Ustop 

CALL TIMCPU(ICP) 

ICP(2) = ICP(2) - ICP(l) 

IF(ICP(2)-5000 .LT. 0) IYES=1 
Subroutine Uthird 

C0MM0N/JP0SMS/FILL(36),ZS(4),ZL(4) 

COMMON/CPOEMS/XMU , WMU , WMU I N , XMUM , XMUS 
Subroutine Uxtoa 

COMMON/T I NDPA (DUMMY 1 ( 4 ) , X I NCRE , DUMMY2 , NDOES 
DIMENSION Y(124) 

Subroutine Vecele 

COMMON/CPEMS/XMU , WMU , WMU I N , XMUM , XMUS 



APPENDIX C. ADDITIONS TO SIMULATION 



Main: Integration additions 

COMMON/GRATON/ITEG 
Subroutine Rbegin 

COMMON/GRATON/ITEG 
READ(5,*)ITEG 
PRINT *,ITEG 
Subroutine RK78 

COMMON/GRATON/ITEG 

IF(ITEG ,EQ. 1) THEN (use variable step size scheme) 

ELSE 

0T6 = DT/6.D0 
DT2 = DT/2.D0 
CALL UDER(T,X,F1) 

T = T + DT2 
DO 200 1=1, N 

DUM (I) = X(I) + DT2*F2(I) 

200 CONTINUE 

CALL UDER(T,DUM,F2) 

DO 210 1=1, N 

DUM(I) = X(I) + DT2*F2(I) 

210 CONTINUE 

CALL UDER(T,DUM,F3) 

T = T + DT2 
DO 220 1=1, N 

DUM(I) = X(I) + DT*F3(I) 

220 CONTINUE 

CALL UDER(T,DUM,F4) 

DO 230 1=1, N 

X(I) = X(I) + DT6*(F1(I) + F4(I) + 2.D0*(F2(I) + F3(I))) 
230 CONTINUE 

ENDIF 

Subroutine Ubegin: 

COMMON/GRATON/ITEG 
READ(5,*)ITEG 
PRINT *,ITEG 

DEPLOYMENT LOGIC 
TL = TLNINI + TLNINl + TLNIN2 

999 IF(TLN .GE. TL) THEN 

IF(IDEPL .GT. 1) THEN 
WRITE(6,*) 'TETHER HAS BEEN DEPLOYED* 

IDEPL = 1 
DVN = O.DO 
DAC = O.DO 
TLNINI = TL 
ENDIF 
TLN = TL 
ENDIF 

IF(TLN .LE. l.D-5) THEN 
IDEPL = 1 

WRITE(6,*) 'TETHER HAS BEEN RETRACTED' 

ENDIF 



APPENDIX D 


TEKTRONICS PLOTTING CODE 
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ill## • FLABY<2#)«FtPt3,Si»FLABXP<8>,FLAiVP(8)# 

• la## f FLAB>CT(8i.FLAiYT(D>»FLAiXlC8)*FLAiyi(8) 

#13## C 

•14## CC ACAD IN PtOTTlNG UAAIABLCS FPOfl THE TCTHCA SINULATION 

•15## C 

#18## OPCN<UfUT-UFlLE -npUil .dai' ^FOPn-^UNFOANATTCD' > STATUS-^ OLD' ) 

tP 

• 17 ## 0PCN(UH1T-2.FILE • ' ip Ui2.d%i' .STATUS-' OLD' ) 

#18B# 0P€N(UNIT-3.FILE -'pploil.dtti' .FOAM-' UrrOAnATTEO' .STATUS- 'OLD' ) 

• 19## 0PEN<UNIT-4.F1LE - *^pp I ol2 .dfti' .STATUS-' OLD' ) 

•2### C 

•81## IP - # 

•82## lU - • 

•23## PEAD12.S) ICOL. IROU 

•24## DO 1# I-l.ICOL 

•as#e RCAD(2. 15) (FLU. J).J-1.8) 

•26## 1# CONTINUE 

•a?## DO 2# I >1. IROU 

•28## READ(l) (X< J.D.J-UZCOU 

#29## 2# CONTirRJE 

#3### READM.Si IC0L2. IR0U2 

#31## 00 21# 1-I.XC0L2 

•32## READ(4.1S) rrLP(I.J).J-t.8i 

*P 

•33## 21# CONTINUE 

•34## DO 22# I-1.IR0U2 

•36## REA0(3)(XP(J.t).J-t.lC0L2) 

#36## 22# CONTINUE 

•37## 1 URITC(S.S) 'UARIA8LC LISTt' 

#38## DO 3# X-l.ICOL 

#39## URlTE(fi.2S) I. (FL( I . J) . J- 1 . 1) 

#48## 3# CONTINUE 

•41## URITE<6.*)'UMAT TVPC OF PLOT#' 

#42## URtTE(6.«>'l-RC0ULAA. 2-UALKlNG. 3-TCNSION. 4-NOK' 

#43## READ( Si 8 )PLTFLQ 

#44## IFIPLTFLQ .EO. 2) GO TO 2#1 

#4### IFCPLTFLG .EQ. 3) 00 TO 3#1 

•44## IFIPLTFLO .EQ. 4) 00 TO t#2 

#47## DO 4# I-liXCOL 

#48## IY(Z) • • 

8 


I 


OF POOR QUALfTY 




C 


•MM 

Itl 

UKm(€#f>'OKC PtaTyM<#H' 

•MM 


UKmcS^tl'CKTCK X-IAKL 19)' 

!•••• 


RCMCS*15mFLMX1< I >,1-U8> 

{•IM 


URm tC#f>^EHTC8 V-UKL <8)' 

tMM 


ReMM5*lS)(FLMVi<I).l*U8> 

1#3M 


CAUL ZHirrcM#} 

1#4M 


CALL TeRnO#1924) 

t#5M 


CALL CRUPLT( 1 « IROU*rLK.FLVt .S^FLAIXI »8*FLA9Y1 ) 

!•••• 

C 

GO TO 1 

!•••• 

I#9M 

1119« 

iiaM 

ii3e^ 

!•• 

URITtf€#«)'TUO PLOTS/CF9F«' 
UIITCCC«t)'eilTCR X-LMCL (29)' 
RCAIMS#35MFLA9X(Z>»1-1#M) 
URmcC^O'CKTCR V-LAlCL '29)' 
REM(8*35MFLA8Y(n«l*UM) 


CALL lNirrr969) 



CALL TCPtm.1924) 



CALL CR^2< 1 « lROU*FLX.rLVl ,FLV2,29.FLA8X*28«FLAr/) 

it6^# 

117M 

C 

GO TO 1 

118M 

193 

Uf»ITEf6«>)' THREE PLOTS/GRAl>H' 

ime 


URITE(6*>>'CMTEP X-LAIEL <29)' 

18M8 


READr5«3Sl(FLRBXCl )»Z-UM) 

121M 


URITErG.O'CHTCR V-LABCL (29r 

128M 


RCAD(S.3S) (FLABYCXiW-IfM; 

12388 


CALL ZNITT<9€9» 

I248# 


CALL TERn(3*l924) 



CALL CPU3(l,lR0U«FLX«FLVl*FLV2.FLY3»29»FLA8X#29»rLRIV> 

126M 

127M 

C 

GO TO 1 

128M 

tP 

I29M 

194 

UPtTE(C«t)'FOUR PLOTS/GRAPH^ 


URITECC.tl 'ENTER X-LA9CL <29)' 

138M 


READ(5*3Sl(FLA|X<n«l-l#B9> 

131M 


URZTECC/O' ENTER V-LAIEL ^29)' 

132M 


READ(S^36) (FLAiV(X)«!-l*2«) 

133M 


CALL XN1TT(969) 

134M 


CALL TCRriO«1924) 

135M 


CALL CAV4a»lR0U,FLX.FLVl,rLV2,FLV3,FLY4,29#FLA»X#89,FLA8V) 

I38M 

137M 

C 

GO TO 1 

138M 

CCC 

UAUCXHG POSITION PLOTS 

139M 

C 


IMM 

Ml 

CONTINUE 

I41M 


IRRU -XROU 

tMM 

U1 

URXTC(8«8>'CHTER NUNRCR OF TETHER NODES' 



RCADC8« 8 INODES 



Hi • NODES ♦ 1 


s 


ORIGINAL PAGE iS 
OF POOR QUALITY 





\ 


c . 'i 

T Xi 


TJ Q 
O 2S 


O 'X* 


50 r- 


lO "XI 




! 








p 


i9im 

C 

ZMX - XP2C3*m#l-l> -- XPCOrNitn 

1P409 

C 

2fm • #i#s(znx) 

UM% 

c 

XAl - ##S<Kn)0 

1W#« 

c 

VA# - A#5(VlfX) 

197«« 

c 

irocpx .or. Atsix^oi xx - xnx 

lOSM 

c 

IFCXnX .LT. #•#) THCH 

19990 

90000 

c 

c 

1P(XAB .OT. AtS(XX>) XX - SIGNfXAi.XnX) 
CHOIF 

90190 

c 

ircvrtx .CT, Aisrvx)) vx • vfix 

c 

IF<Vnx «LT. #•#> THCN 

90300 

c 

inVAi .CT. AfS(VX)) vx • 51QN(VAi*Vf#0 

90400 

c 

CNDIF 

m^^WW 

c 

iFunx .or, Aisczxi) zx • znx 


c 

IFiZnX .LT. #«#) THCN 

t«7#« 

c 

IFrZAi AtSIZX)) ZX • SIGN! ZAt«ZAX) 

90000 

tp 

c 

CHfilF 

2e9e« 

C 4## 

CONTll«je 


c 

URITC(6*S>'XX*VX*ZXf ^«KX»VK#ZX 

21122 

c 

ZX2 - .S«2X 

91900 

c 

1F(YX .CO. #.#> GO TO 4#1 

91300 

c 

VSCAL - AlSfZXZ^VX) 

91400 

c 

00 TO 4*2 

91000 

C 4#| 

VSCAL » 1-# 

91000 

C 4#2 

COMTINUC 

9X700 

c 

XSCAL - AB$(Zxa/XX) 

sisee 

cc 


9X900 

cc 

FIND HAXinun DCFLCCTION IN THE X AND V AXES 

99000 

98X00 

cc 

XF(1U .EO. #> THEN 

99900 


XX • #.# 

99300 


VX • #.s 

99400 


HINT • 0 

tP 

99009 


DO 4## I-l.tARU 

99000 


DO 4#2 J-2.NOKS 

98700 


XNX • XPU,J4NINT) - XPa.I^HlNT) 

99900 


VHX • XP(8*J4N1NT) - XP(2*1^N1NT) 

99000 


XAtS - AB9(xnx> 

93000 


YABS - AiS(VflK) 

93X00 


1F(XA»S ,Gt. XX) THEN 

93900 


XX • XAM 

93300 


IXX • I 

93400 


JX • J 

93000 


xstoN • xm 

93009 


ENDIF 

93700 


IFtVAB# .OT. VX) THEN 

03000 


VX • VAii 

93900 


IVY • X 

24### 


JV • J 


Onf^NAL PAGE !S 
OF POOR QUALITY 



vnx 




f41#« 



YflON • vnx 

|4Mt 



CHOir 


408 


corrriNuc 




«IMT • NIKT ♦ Nt 


400 


COtiTINUC 

E4C«« 



U81TC<0*t>^n0X IN PLANE DCFUECTION • ' »X5XGI« 

«47#0 



URlTC(0«trAT THE' »1XX.' TINE STEPj NOC€ 8' * JX 

94WNP 



iMITE(6«tri1AX OUT or FLAHE OCFUCTION -',VS1GN 

t490# 



URITE<0«*)^AT THE',1W,'TINE STCPj NODE 8' , JV 




lU • 1 




CNOir 




N48 • 0 

«S3ee 



N38 • 1 




U81TECC»8>'ENTER XSCAL AND VSCAL' 

tSSM 



R£A0( $ « 8 ) XSCAL » VSC AL 


C 


U8I TE ( 0 »«)' XSCAL » VSCAL 1 ' »XS€AL* VSCAL 





2S7M 



DO 430 l-l,XffOU 




DO 440 J«1,N0DE5-1 




JJ • J^1^N48 


c 


XP2(l«JJi - XSCAU<XP8(UJJ) - XP2(1*NU1)) ^ 

SilM 

c 

8 

XP2(l.NltI) 

2C2O0 

c 


XP2(2,JJ) • VSCALKXP2(2,JJ) - XP2(2.N18in 4 

a63«« 

c 

8 

XP2(8»N18I) 

a«4#e 



XP2<1#JJJ - XF2(1*JJ) - XP2(l,NttI) 4 XSCAL8XP2(1 

26SM 



XP2f2«JJ) • XP2(2*JJ> “ XP2(2«Hltl) 4 VSCALtXP8(2 

244## 

440 

CONTIflUC 

247M 

C 


URirE(6»8)'XP2(l*3)#XP<1.3)l' »XP2( i *3) »XP( 1 « 3> 

e<8M 



N42 • N42^m 

B6944 



XP8(lWi32> • X$CAL8XP2(1#N32> 

87404 



XP2( 1.N324N0DES> • XSCALtXP2(l*N324N0DCS> 

87100 



XP2(2,h32) • VSCAL<XP2(2.K32> 

87800 



XP2<2«N324H0DES) • YSCAL8XP2(2«N324N0DES) 

87300 



H32 - N32 4 N1 

87400 

430 


CONTINUE 

87S00 

C 


endif 

87400 



UA1TE(6*«I'UARIA8LESi' 

87700 



U81TE(0»*I MN-PLAHE, OUT-OT-PLANE • RADIAL' 

17800 



UR1TE(0#K) 'ONE PL0T/G8APM' 

87900 



URITE(0*X) 'INPUT X AND V UARIABLE NUN8EII' 

88000 



RCADCSrt) IXP«IYP 

88100 



DO 838 I-1«IC0L8 

88800 



1F<1XP .CO. I) THEN 

88300 



DO 838 J-I.IR0U8 

88400 



FLXP(J) • XP8<I#J) 

88S00 

838 


CONTINUE 

88800 



DO 230 K*l#8 

88700 



FLAtXPCK) - rtJP(I*K) 

88S00 

838 


CONTINUE 


« 


OF POOR QUALrnr 






Mi00 

3112#^ 


393«# 


fP 

305M 

M9«# 

31m 
3I1M 
312M 
313M 
31 m 
315M 
3169# 
31709 
31800 
31900 
32000 
tP 

32100 

32200 

32300 

33400 

32800 

32800 

32700 

32800 

38900 

33000 

33100 

33200 

33300 

33400 

33900 

33000 

8 


ccc 

c 

301 


310 


311 


330 

320 


300 

340 


202 

C 

f 


CNDIP 

iruvf .CO. 11 THCM 
DO 240 J*t*10OU8 

FLYP<J) • XP2<1«JJ 
COffTXHUC 
DO 241 K»l,n 
FLABVPfCl - FLP<l»ICi 
COHTIHUC 
CHDIF 
CONTXNUC 
C0LL IfllTT<900) 

CALL TCPTfO, 10241 

CALL CAUPri»IfK)y2*FLXF.rLVP,8»FLABXF.8»FLAtVP*N0DCS#IIKlU 
IPOU • iimi 
GO TO 1 


TCNS10T4 PLOTTING CODC 
COMTINUC 

URITC<6.jn 'ENTER 8ASE NUN8ER OF UARXAIIXS' 
READ(5#>)I8ASC 
HAINl • ICOL - I EASE 
DO 310 J-WIROU 

FLKTrj) - X<1»J) 

CONTXNUC 
DO 311 K-1.8 
FLAiXT(K) • FL(t*IC) 
rLAIVT<K) - FL<18ASC«1.IC) 

COKTINUC 

DO 320 1- IBASC4U1C0L 
II • I - IlASe 
DO 330 J«1*1R0U 

VOATdl.Jl • X(UJ) 

CONTINUE 
CONTINUE 
DO 340 1-l.NNZNl 
II - I-l 
DO 3S0 J'l.IROU 

V(J^1I81R0U1 • VDAT(I.J) 

CONTINUE 
CONTINUE 
CALL IN1TT(O00> 

CALL TERm3d024> 

CALL TCHCRO<l,IROU.NRlHI»FL>a,V,8,FLA8XT,8#FLArm 
GO TO 1 

12I1TC(0«8) 'FINX8HCD PLOTTING' 

FORNATCtlfl 








OF POOR QUALrTY 






^ -lb -5 


9010m 

C 



00000 

ccc 


SUiPOUTINC TO PLOT 2 CURUCS ON 09C CRPPH 

00300 

00^00 

c 


SU920UT1NC CmAe2CIH*lF0»XDPTP,YIMTAI»VlMTA2.LENX, 

00000 

c 

0 

LEHV.LAIV) 

00700 



DINENSIOH Xl>ATA<2O#O>»VDATAt(2t00)«yiMTA2fa#«0l* 

00000 

00900 

c 

9 

LABX(2#)«LAlVCat) 

01000 



ifllN • X0000. 

01100 



iHAX • -X0000. 

01000 



M1IN • BIUN 

0X300 



ANAX • iflAX 

0X400 



CALL ilNITT 

0X000 



CALL CHRSlZfZ) 

0X000 



CALL NPTS(IFO) 

kP 

0X700 



CALL XFPf1(2) 

0X000 

0X900 

c 


CALL VFRA(2) 

00000 

ccc 


FIND niNlNUn AHD NAXinUR UALUCS 

eaiM 

00B00 

c 


CALL nNnX<YDATAl«tnjM«inAX> 

00300 



CALL hNnX<VDATA2.if1]N*BnAX> 

0Z400 



CALL nriMXCXDATA.AniH/ANAX) 

00S00 



URITC(6.<rxPIN«' .Af1|N,'XnAX-'«AI1AX 

00600 



URlTC(6.f )'UAKT NCU X-AXIS LIMITS71-VtS»2-N0' 

03700 



READCS.kULiX 

03900 



IFriLXX .NE. U GO TO 1$ 

0Z900 



URITE(6«t)^CNTER xmN#Xf1AX' 

03000 



READ<S»«)AniH*ARAX 

02X00 

15 


CALL DLinX(AI1tN*ArfAX) 

02300 



URITE(6**)'VniW',i«IH#'V«AX*',W1AX 

IP 

03200 



URITCr6,»)'D0 VOU UANT NCU V LlfllTST 1-VtS# 2-MO 

02400 



READ<5«f )ILin 

02S00 



IFULIM .NC. 1) GO TO !• 

02600 



URITCCS^O' ENTER VflXN AMD VflAX' 

02700 



RCAD<5«t)BniN»DMAX 

02000 

15 


CALL DLXMVfBNXN.iNAX) 

02900 



CALL CHR5lZO> 

04000 



CALL CHCCIC(XDATA.VDATA2) 

04X00 



CALL DSPLAV<XDATA«YOATAa> 

04300 



XX - <750 - LCNXS13)/8 « 155 

04200 



IV • <575 - LEffVkill/t 4 125 ♦ LCNVt2l 

04400 



CALL CHR$IZ<gl 

04000 



CALL n0UAi5<XX.25) 

04000 



CALL TSCMD 

04700 



CALL MLAiCLCLEflX^LAtX) 

04000 



CALL 5XZE5(.£SI 


t 


OF POOR QUALITY 



wuw?!v!IfIiIiuIu55i!*-*^*-SSIlSTlS*-S5#S5555««®888** 

mmiilllHiH iHIUmiilelllS liillllliillillli 


IS 


IS 



SCA0<5.«>rLlX 

IFIILIK .fC. 1> 00 TO IS 

USITCCfi.trCKrCS XffIff,)«l%K' 

SCAO< 5 * t > SniH * AfMX 

CALL DLlKXKAniN.ATMX) 

usiTC(S^*> 'vmN-',sni«,'VfiAx-',inAx 

mUTCfS*t)'E>0 vou UAHT HCU LIKITS ON THE V-AXIS? 1-VCS»2-N0' 

REAOCS*S) ILIN 

IFCrLin .1C. 1/ GO TO If 

UffirC<6*ti'EHTCS V-^NIN AHO V-f1AX' 

FCA0<5»f) fniN«BfSOC 
CALL DLiny(tNXN«tt 1 AXf 
CALL CHPSIZO) 

CALL CHECKfyDATA«YDATA3) 

CALL 0SSLAV(Xt>ATA.Vt)ATA3f 
IX * 4[7SS - LCNX413//2 4 ISS 


IV • '57S - LEMV321I/2 ♦ 125 ♦ LENY<21 
call CHSSIZ(2) 

CALL »0MAfS(IXi>25J 
CALL TSCND 

CALL HLABCLCLCNX.LABX) 

CALL SIZES( .25> 

CALL SVABL(120> 

CALL STEPSCIS) 

CALL CPLOT(XOArA«VOATAl> 

CALL S'/nBL(l21> 

CALL STEPSdS) 

CALL CPL0T(XDATA»VDATA2> 

CALL CHRS1Z<2I 

C CALL nOUABS(2SS*72S> 

C CALL HLABtLCLENV^LABY) 

CALL nOUABS<25»IV) 

CALL ULABEL(LENV.LASY) 

CALL CMASIZ(>n 
GO TO U«2) IH 

2 CALL riNlTT(S#7SS) 

CALL HONE 
READIS^T) 

call NEUPAQ 
GO TO 3 

1 CALL MDCOPV 

CALL NEUPAG 
CALL FINITT(S«7SS> 

3 RETURN 
END 

C 

ccc 

CCC SUBROUTINE TO PLOT 4 CURUEf ON ONE GRAPH 


ORIGINAL PAGE IS 
OF POOR QUALITY 



II 


• 84 M 




• 62 M 


• 64 «« 

tP 

665 M 


•C» 80 « 

• 89 M 

•7000 

•7100 

0780« 

07300 

07400 

07S00 

07000 

07700 

07800 

07000 

00000 

*P 

00100 

00800 

00300 

00400 

00800 

00000 

00700 

08000 

00800 

•9000 

00100 

00800 

•0300 

09400 

09000 


C 

C 


8 


i 


3 


C 

ccc 

ccc 

ccc 

c 


c 


c 


c 

c 

c 


CALL fvraui20) 

C0LL 0T8FS(101 

CALL CP1OT(X©0TA,VI>0701) 

CALL CMRSIZISJ 

CALL nOOAiS<i^.>72S> 

CALL HUA0€L(LCrfV«LA8Y) 

CALL rNMM0S<2S.IY) 

CALL MLA0€L<L£NY*LA0Y) 

CALL CH0S1Z(4I 

GO TO <1#8>1H 

CALL riHITTC 0,7001 

CALL Horc 

R£AD(S,n 

CALL NCUPAG 

GO TO 3 

CALL H0COPY 

CALL fCUPAG 

CALL F2N1TT( 0,7001 

RCTUPH 

EKD 


SUlPOOrilC TO PLOT 3 CURLCS ON ONE ORI^ 


SUBROUTINE CPU3< IH. IFO/XMTA, VDATAl , YDATA2,VDATA3,LENX, 

0 LA8X,LENY,LABV^ 

DinCNS I ON XDATA( 8000 ) , VDATA 1 ( 2000 > , YDATA2 C 2000 ) , VDAT A3 < 8000 ) , 
8 LABX(80>,LAIV(201 

iniN - 10000. 

BNAX • - 10000 . 

ANIN - BfUN 
ANAX • BNAX 
CALL 8INITT 
CALL CHRStZcai 
CALL NPTS(IFO) 

CALL XFRn(2) 

CALL vrRN(2) 

FIND niNinun and naxinun ualucs 

CALL nNnx(yDArAi,0mN,inAx) 

CALL f1NnXiyf>ArA8,BN2N,BNAX> 

CALL nNnX(YDATA3,8NlN,>NAX) 

CALL nNNX(XDATA,ANIN,ANAXl 

URXTC(0,t)'UANr NEU X-AXIS LIHITSTl-VES.R-NO' 




OF POOR QUALITY 




149## 

IS100 

1S3«# 

1S4M 

19990 

1S6M 

1$7«« 

ism 

1S9M 

16999 

*P 

UlM 

1624# 

16340 

16440 

16540 

16644 

16740 

16644 

10944 

17440 

17140 

17240 

17344 

17444 

17544 

17644 

»P 

17744 

17444 

17940 

14440 

14144 

14244 

14344 

14444 

18540 

14644 

14740 

14440 

14940 

19440 

19140 

19000 


\ 


ccc 

c 

4U4ITOUTINC CWV41 ]H. IPO.KMTOr V9PTOW VlMT02«yiMTA3, VOPTM* 

4 UCr«<,LAIX*IXr<V»LA4Vl 

c 

Dlr«f«S10H XDAT4(2#00)*VDATA1I2000>*VDOTA2(2000I»VDAT43(2000)* 
4 VlMT04<2000>*L4»X<20)*LA4Vf20> 

C 

4ntN • 10440. 

OfVkX • -14400. 
oniN - 4P1N 
onox • onox 
COtt 41NITT 
CALL CHPS1Z(2) 

CALL HPTSaFOl 
CALL XFR11(2) 

CALL VF«n( 2 ) 

C 

C FIND niNinUPI AND NAXinutf UALUCS 

C 

CALL nNNX<VDATAl»4firN»t^> 

CALL nNNX(VDATA2.BniN*illAX) 

CALL NHNX(VDATA3«BmN.BnAXi 
CALL m<nX(VDATA4»INlN«inAX) 

CALL nMnX(XDATA*AHIN/AnAX) 

URITC(6.X)'Xf1IN-' .ANIM^^XNAX*' .AhAX 
URITE(6*t)'UANT N£U X-AXI5 LimTS71-ytS#2-H0' 

PEAOrs^OUlX 

IFIILIX .NE. 1) GO TO IS 

URIT€(6#f >*ENTCR XHlN^XfMX' 

ReAD(5«t)AniN.ANAX 
IS CONTINUE 

URITE(6*«) 'BNIN-' ^OniN.^BNAX** .EfIAX 

WRlTt(6#»)'D0 YOU UANT NEU LINITS ON TWC Y-AXIS? l-VES,2-«0' 

READ(S*B) ILIfl 

IF<lLin .NE. n GO TO 14 

UR1TE(6«S)'ENTCR V-fllH AND V-tlAX' 

RCAD($.t) BRIN^BNAK 
10 CALL DLlNY<8nlN.BnAX) 

CALL CHRSIZO) 

CALL CHECK <XDATA«V|MTM} 

CALL DSPLAV<XDATA«VOATAU 
IX - (754 - LENXB13I/8 ♦ 160 
lY • (STS - LCNY48n/8 « 12S ♦ LCNY421 
CALL CH«^1Z<2) 

CALL N0UA4S(1X*06) 

CALL T9CN0 

CALL HLA8tL(LSHX*LAiX> 


te-»a 

L'# 
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OF POOR QUAU" 
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OF POOR 


quality 


I 


3 

F - 

i 

> 


S 

•- 

s 


k 
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DO f l«L*Mlim 
II • I -1 
DO IS J-I.IFO 

Tcnpcj) - v<j4ii»iro> 

COHTlHUe 

CALL nNrtx(rcnp«BniN«inM() 

CONTINUE 

CALL nNAX<XDATA»AmN*AnAX> 

UPITE(6*BI'XmH-' ,«1IH,'X«AX-',ANAK 
URITCf€.t>'UAKT rCU X-AKIS LIHITS7l-VCS#2-N0' 

AEADI5#f )ILtX 

IFdLlX .NE» 1) CiO TO 8 

UffITCI£«t) 'ENTER XHIN^XNAX' 

READ< 5 » « ) ANtN* ANAX 

CALL DLINX^ANlr^AflAX) 

URlTE(6*tl 'VniN-'.iNlH.'YnAX-'^BMAX 

URITE(6«t)'tK) VOO UAHT NEU LINITS ON THE V-AXIS7 1-YES#2-^ 

READ(S«tl ILin 

irciLin .NE* i) GO TO 18 

UR1TE(6*S)' ENTER V-RXN AND V-IMX' 

READ<5#8> BniN.BNAX 
CALL DLinY(BniN,BnAX) 

CALL CHRSIZOI 
DO 2S J-l.IFO 

TENPU) • V( j4lFO«<NniNl-m 
CONTINUE 

CALL CHECK (XDaTA.TENR) 

CALL DSPLAV<XDATA,TEnP) 

IX • (758 - LENX413)^2 ♦ 158 

IV • (57S - LENyt21l/2 ♦ 125 4 LENVB21 

CALL CHRSIZ<2) 

CALL notJABS(lX*2S) 

CALL TSEMO 

CALL HLABEL(LENX#LABX> 

CALL SIZES(.26) 

DO 35 I-i«NNINl-l 
II • I-l 
DO 85 J-ldFO 

TENP(J) - Y(j4lXtlF0) 

CONTINUE 

CALL CnOT<XDATA.TEnP) 

CALL SVnBLC 11941) 

CALL STEP9(54BXI> 

CONTINUE 
CALL CHRSIZIB) 

CALL N0UAB8(M8*78S> 





fp 

306e« 

3«7M 

3 « 80 « 

3i«e« 
3iie« 
31 200 
31304 
31400 
31000 
31600 
31700 
31800 
31900 
32000 
»P 


32100 

32200 

38300 

32400 

32600 

32600 

38700 

38100 

38000 

33000 

33100 

33800 

33300 

33400 


33000 


33600 


8 



C CALL HLAECLOXNV.LAiV) 

CALL nOUAR8(2S*lV> 
call VLA86LUXHV*LA8Y> 

CALL CHR01Z(4> 

GO TO (I.2J IH 
8 CALL riNITT(0,7001 

CALL HOfie 
8CA0(5.8> 

call NCUPAQ 

GO TO 3 

1 CALL HOCOP*/ 

CALL HCUPAG 
CALL FtNtTT(0«700) 

3 RETUPH 

CND 

SU8POUT1NE CR^0’<lH#NPTS9,XARlMV»YARRAV«LCH><»LA8X«LErfV*LA8V<i 

8 NODES. IROU) 

C 

C SUDPCUTIfC TO PLOT UAUCING CURUCS OH A GRAPH 

C 

DinCNSlOH XARPAV(7000).VARRAY(7000J#LA8X(8I#LA8V(8) 
DlfCHSIOH XTtf»*l7000i.yTE«P(7000) 

C 

AHIN • 10080. 

AHAX - -10000. 
iniN • AHIH 
8HAX • AHAX 
CALL ilNITT 
C 

call CHR51Z(2) 

CALL NPTS(NPTSSl 
CALL XFRm2> 

call vri?fi<2) 

CALL PN«X<XARRAV.A«IH.AhAXl 
CALL nNNX(VARRAY.iHlN.Bf0OO 
W«ITi:(6,»)'XF1lN-'.AHIN.'X«AX»^ .AHAX 
URITt(6.t)'UAKT NCU X-AXIS L1MT$7 1-VCS.2-HO' 

RCAD(5«f >1LIX 

IFCILIX .NC. 1) 00 TO IS 

URITCCS.Sl'CNTER XniN.XfMX' 

READ(S.81AmH.AHAX 
iS call 0LIPX(AHIN.AHAX> 

URI TE ( 0 . 8 ) ' vn 1 H - ' . ilU N « ' VnAX - ' , iflAX 
URITE(«.8)'UAHT NEU Y-AXIS LIHITS? 1-VES.8-N0" 

READ(S.8)ILIV 

IF<XL1Y .HE. n GO TO 80 

UR1TC(0»8)'CNTER VHlH.VfVIX' 

REA0(f*8)8mN.8tMX 
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